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1 .  INTRODUCTION 


Aerosols  commonly  found  at  stratospheric  altitudes, 
10  -  30  km,  are  a  result  of  photochemical  formation 
involving  sulfur  compounds.  These  background  aerosols  are 
generally  uniform  over  the  globe.  The  concentrations  of 
these  stratospheric  aerosols  can  be  increased  dramatically 
following  massive  volcanic  eruptions.  Additional  sulfur 
based  aerosols  can  be  photochemically  created  from  the  sul¬ 
fur  gases  in  an  eruption  cloud.  Volcanic  ash  and  debris  can 
also  be  injected  into  these  altitudes. 

The  Air  Force  Geophysics  Laboratory  (AFGL)  has  devel¬ 
oped  three  stratospheric  aerosol  models^  for  use  in  the 

2  3 

transmittance/radiance  models  LOWTRAN  and  FASCODE  .  The 
three  models  are  for  background  stratospheric  conditions, 
fresh  volcanic  aerosols,  and  aged  volcanic  aerosols.  These 
models  have  also  been  adapted  for  use  in  the  Standard 

4 

Radiation  Atmosphere  (SRA)  models. 


1.  Shettle,  E.  P.  and  R,  W.  Fenn  (1976)  "Models  of  the 

Atmospheric  Aerosols  and  Their  Optical  Properties,"  in 
AGARD  Conference  Proceedings  No.  183,  Optical 
Propagation  in  the  Atmosphere.  Presented  at  the 
Electromagnetic  Wave  Propagation  Panel  Symposium, 
Lyngby,  Denmark,  27-31  October  1975,  AGARD-CP-183 , 
ADA028-615. 

2.  Kneizys,  F.  X.,  Shettle,  E.  P.  ,  Gallery,  W.  O., 

Chetwynd,  J.  H.,  Abreu,  L.  W. ,  Selby,  J.  E.  A., 

Clough,  S.  A.,  and  Fenn,  R.  W.  (1983)  Atmospheric 
transmittance/radiance:  Computer  code  LOWTRAN  6, 

AFGL-TR-83-0187 ,  AD  137796. 

3.  Clough,  S.  A.,  Kneizys,  F.  X.,  Shettle,  E.  P. ,  and 

Anderson,  G.  P.  (1986)  Atmospheric  radiance  and 
transmittance:  FASCOD2 ,  Proceedings  of  the  Sixth 

Conference  on  Atmospheric  Radiation,  Williamsburg, 
Virginia,  American  Meteorological  Society,  Boston, 
Massachusetts,  141-144. 

4.  WCP  (1986)  A  preliminary  cloudless  standard  atmosphere 

for  radiation  computation,  WCP-112,  WMO/TD-No  24, 
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The  background  stratospheric  aerosol  consists  of  a  75  % 
solution  of  sulfuric  acid  in  water.  The  wavelength  depen¬ 
dent  index  of  refraction  for  the  solution  is  based  on  labo- 

C  £  7 

ratory  measurements  at  300  K  '  '  .  The  indices  of  refrac¬ 

tion  for  the  volcanic  ash  are  based  on  measurements  of 

Q 

Volz  .  The  size  distributions  for  all  three  stratospheric 
aerosols  are  given  by  modified  gamma  distributions. 

Stratospheric  temperatures  are  well  below  300  K.  Val¬ 
ues  of  the  index  of  refraction  for  H_SO.  are  now  available 

2  4 

o 

at  250  K  and  can  be  used  to  calculate  temperature  dependent 
indices  of  refraction  for  for  temperatures  appropriate 

for  stratospheric  conditions.  Also,  recent  measurements  of 
the  size  distribution  of  stratospheric  aerosols  indicate 
that  a  log  normal  distribution  may  be  more  appropriate  for 
stratospheric  aerosols  than  a  modified  gamma  distribution10. 


World  Climate  Research  Programme,  International 
Association  of  Meteorology  and  Atmospheric 
Physics/Radiation  Commission. 

5.  Remsberg,  E.  E.  (1971)  Radiative  properties  of  several 

probable  constituents  of  Atmospheric  Aerosols,  PhD. 
Thesis,  Department  of  Meteorology,  University  of 
Wisconsin,  Madison. 

6.  Remsberg,  E.  E.  (1973)  Stratospheric  aerosol  properties 

and  their  effects  on  infrared  radiation,  J.  Geophys . 
Res. ,  78:1401-1407. 

7.  Palmer,  K.  F.  and  D.  Williams  (1975)  Optical  constants 

of  sulfuric  acid;  Application  to  the  clouds  of  Venus?, 
Applied  Optics,  14:208-219. 

8.  Volz,  F.  E.  (1973)  Infrared  optical  constants  of 

ammonium  sulfate,  Sahara  dust,  volcanic  pumice,  and 
flyash,  Appl .  Opt . ,  12:564-568. 

9.  Pinkley,  L.  W.  and  Williams,  D.  (1976)  The  infrared 

optical  constants  of  sulfuric  acid  at  250  K,  J,  Opt. 
Soc .  Am . ,  66^122-124. 

10.  WCP  (1983)  Report  of  the  Experts  Meeting  on  Aerosols  and 

Their  Climatic  Effects  ,  Eds.  Deepak,  A.  and  Gerber, 

H.  E.,  World  Climate  Programme,  WCP  55,  World 
Meteorological  Organization,  p  30. 
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1 . 1  Organization  of  Report 

This  report  describes  the  development  of  a  new  back¬ 
ground  stratospheric  aerosol  model  in  which  temperature  de¬ 
pendent  indices  of  refraction  for  have  been  used  with 
a  log  normal  size  distribution.  Chapter  2  is  a  summary  of 
the  literature  that  provides  the  justification  for  the  de¬ 
velopment  of  a  new  model.  Chapter  3  describes  the  methodol¬ 
ogy  used  to  formulate  the  new  model.  Chapter  4  presents  re¬ 
sults  from  the  new  formulation  and  compares  them  against  the 
formulation  currently  used.  Finally,  Chapter  5  presents  the 
conclusions  and  recommendations  for  future  work. 
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2.  A  REVIEW  OF  STUDIES  OF  STRATOSPHERIC  AEROSOLS 

Two  recent  major  volcanic  eruptions  have  provided  re¬ 
searchers  with  the  opportunity  to  study  stratospheric  and 
volcanic  aerosols  in  depth.  The  eruptions  were  Mt.  St.  He¬ 
lens  in  1980  and  El  Chichon  in  1982.  In  both  cases,  the 
eruptions  were  intensely  studied  by  ground  and  aircraft 
based  lidars,  balloons,  aircraft  and  satellites. 


2 . 1  Lidar  Studies 

Lidar  provides  a  tool  to  study  particulate  matter  in 
the  stratosphere.  Lidars  at  various  locations  around  the 
world  have  permitted  researchers  to  study  the  distribution 
of  volcanic  materials  injected  into  the  stratosphere  since 
the  1970's.  The  altitude  resolution  that  can  potentially  be 
obtained  from  lidars  allows  researchers  to  study  both  the 
temporal  and  spatial  details  of  volcanic  eruption  clouds. 
The  worldwide  network  of  lidars  ( e .g # H > 13 , 14 , 15 , 16 ^  pro_ 


11.  Reiter,  R.,  Jager,  H.,  Carmuth,  W. ,  and  Funk,  W.  (1979) 

The  stratospheric  aerosol  layer  observed  by  lidar 
since  October  1976.  A  contribution  to  the  problem  of 
hemispheric  climate,  Arch.  Met.  Geoph.  Biokl.,  Ser.  B, 

27:121-149. 

12.  D'Altorio,  A.,  Viconti,  G.,  and  Fiocco,  G.  (1981)  Lidar 

detection  of  volcanic  aerosols  in  the  atmosphere 
following  the  Mount  St.  Helen  eruption,  Geophys.  Res. 

Let .  ,  8^:63-65. 

13.  Clemesha,  B.  R.  and  Simmonich,  D.  M.  (1978)  Strato¬ 

spheric  dust  measurements,  1970-1977,  J .  Geophys . 

Res. ,  83:2403-2408. 

14.  Iwasaka,  Y. ,  Hayashida,  S.,  and  Ono,  A.  (1983)  Increas¬ 

ing  backscattered  light  from  the  stratospheric  aerosol 

layer  radar  after  Mt.  El  Chichon  eruption,  laser  radar  * 

Measurement  at  Nagoya  (35°N,  137°E) ,  Geophys.  Res. 

Let . ,  6:440-442. 

15.  Hirono,  M.,  Fujiwara,  M.,  Shibata,  T.,  and  Kugumiya,  N. 

(1981)  Lidar  observations  of  Mt .  St.  Helens  in  May 
1980,  Geophys.  Res.  Let.,  9:1019-1022. 
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vided  extensive  coverage  of  the  eruption  clouds  from  the 


Mount  St.  Helens  and  El  Chichon  eruptions.  These  lidars 
have  also  detected  the  presense  of  so-called  "mystery 
clouds",  clouds  that  cannot  be  traced  to  ■'•he  eruption  of 
known  volcanoes . 

The  lidars  used  for  probing  the  upper  troposphere  and 
stratosphere  use  the  lidar  backscattering  ratio  to  identify 
layers  of  non-molecular  scattering.  The  lidar  back-scatter¬ 
ing  ratio  as  a  function  of  altitude,  B(z),  is  given  as 

B ( z )  =  1  +  fA(z)/fM(z)  (1) 

where  fA(z)  is  the  aerosol  backscattering  function  and  fM(z) 
is  the  molecular  backscattering  function.  A  requirement  for 
the  use  of  (1)  is  the  assumption  that  at  some  altitude  the 
returned  lidar  signal  is  only  from  molecular  scatterers. 
This  altitude,  often  called  the  matching  altitude,  is  used 
to  determine  the  atmospheric  density  profile  required  to 
give  f M ( z ) .  Once  this  is  known,  any  return  greater  than 
that  produced  by  a  pure  molecular  atmosphere  is  assumed  to 
be  from  particulates.  The  matching  altitude  is  generally 
taken  to  fall  in  the  30  km  range,  although  Clemesha  and  Si- 
monich^3  feel  that  the  matching  altitude  should  be  taken  at 
higher  altitudes. 


16.  McCormick,  M.  P.,  Swissler,  T.  J.  (1983)  Stratospheric 
aerosol  mass  and  latitudinal  distribution  of  the  El 
Chicon  eruption  cloud  for  October  1982,  Geophys.  Res. 
Let. ,  9:877-880. 


Table  1  summarizes  results 


12,13,14,15,17,18,19,20,21 


of 


lidar  studies  made  under  background  and  volcanic  conditions. 
The  studies  indicate  that  a  backscatter  ratio  of  between  1.1 
and  1.4  is  representative  of  background  stratospheric  condi¬ 
tions.  The  backscatter  ratio  following  an  eruption  is 
highly  variable  and  depends  on  the  force  of  the  eruption,  as 
measured  by  the  amount  of  material  injected  and  the  height 
of  the  eruption  cloud,  as  well  as  the  sulfur  content  of  the 
eruption  cloud. 

The  lidar  studies  have  detailed  the  layering  of  mate¬ 
rial  following  volcanic  eruptions.  The  layers  are  highly 
dynamic  and  can  change  somewhat  quickly  as  a  result  of  vari¬ 
ations  in  the  stratospheric  circulation  and  as  particles 
settle  out. 


17.  Clemesha,  B.  R.  and  Simmonich,  D.  M.  (1983)  Lidar 

observations  of  the  El  Chicon  dust  cloud  at  23°S, 
Geophys.  Res,  Let.,  10  :  321-324 . 

18  wirono,  M.,  Fujiwara,  M.,  and  Shibata,  T.  (1981)  Lidar 
observation  of  sudden  increases  of  aerosols  in  the 
stratosphere  caused  by  volcanic  injections.  I. 
Soufriere  1979  event,  J.  Atmos.  Terr.  Phy,  43:1127- 
1131. 

19.  Hirono,  M.,  Fujiwara,  M.,and  Shibata,  T.  (1982)  Lidar 

observation  of  sudden  increases  of  aerosols  in  the 
stratosphere  caused  by  volcanic  injections.  II.  Sierra 
Negra  event,  J.  Atmos.  Terr.  Phy,  44 : 811-818 . 

20.  Hirono,  M.  and  Shibata,  T.  (1983)  Enormous  increase  of 

stratospheric  aerosols  over  Fukuoro  due  to  volcanic 
eruption  of  El  Chicon  in  1982,  Geophys.  Res.  Let.,  10 ; 
152-154. 

21.  D'Altorio,  A.  and  Visconti,  G.  (1983)  Lidar  observations 

of  dust  layers'  transience  in  the  stratosphere 
following  the  El  Chichon  volcanic  eruption,  Geophys . 
Res.  Let . ,  10:27-30. 


Table  1.  Summary  of  Lidar  Studies  of  Background  and 
Volcanic  Aerosol  Conditions  in  the  Stratosphere 


LOCATION 

WAVELENGTH 

TIME  PERIOD  ( urn) 

ALTITUDE  REGION 
Uun) 

AVERAGE  BACKSCATTER 
RATIO 

Brazil13'17 
(23°S,  46°W) 

July-August  1970 

0.589 

20 

1.4 

1973 

0.589 

20 

1  06  -  1.12 

August  1975 

0.589 

20 

1.28 

April  1976 

0.589 

20 

1.15 

1982 

0.589 

15  -  20 

>  5.0  (Peak) 

May  1978-Apr  1979 

0.694 

16 

1.052  +  0.02 

May-June  1979 

0.694 

16 

1.12  -  1.42 

Japan 

May  1978-Apr  1979 

0.694 

20 

1.116  +  0.05 

( 33°N,  1 3  0°E) 

Nov  1979 

0.532 

21 

1.05 

(Ref  14,15, 

Nov  1979 

1.064 

21 

1.46  +  0.19 

18,19,20) 

Dec  1979 

0.532 

21 

1.2  (Peak) 

Dec  1979 

1.064 

21 

2.0  (Peak) 

Apr  1982 

0.532 

24.5 

400  (Peak) 

,  20 

Japan 

( 35°N ,  1 3  7  °E ) 

Apr-July  1982 

0.694 

24  -  26 

44  (Peak) 

ICgly12  o 

Sept-Dec  1979 

0.589 

14  -  20 

<  1.2 

(42°N,  23°E) 

June-July  1980 

0.589 

14  -  20 

2.0 

Mar-Apr  1982 

0.589 

15  -  25 

1.5  ( Peak) 

2.2  Balloon  Measurements 

Hofmann  and  Rosen  have  made  extensive  measurements  of 

22  23  24  25 

background  stratospheric  and  volcanic  aerosols  '  '  us¬ 

ing  balloon-borne  instruments  that  can  reach  altitudes  up  to 
about  35  km.  Their  typical  balloon  package  carries  optical 
particle  counters  to  measure  condensation  nuclei  for  parti¬ 
cles  with  radii  r  >_  0.01  ym  and  a  dustsonde  for  particles 

24 

with  r  0.15  and  r  0.25  ym.  A  more  recent  version  also 
carries  a  large  particle  counter  for  r  >  0.25,  0.95,  1.2  and 
1.8  ym. 

The  particle  size  measurements  provide  integral  values 
of  the  particle  concentrations  with  sizes  greater  than  or 
equal  to  the  respective  cut-off  radii.  By  taking  ratios  of 
these  concentrations,  a  set  of  size  ratios  can  be  obtained 
that  can  be  related  to  size  distribution. 

The  balloon  packages  have  also  been  equipped  with  an 
intake  heater  that  heats  the  air  samples  to  around  150°  C. 
This  allows  one  to  distinguish  from  volatile  and  nonvolatile 
aerosols.  This  technique  has  identified  the  volatile 


22.  Rosen,  J.  M,,  Hofmann,  D.  J.,  and  Laby,  J.,  (1975) 

Stratospheric  aerosol  measurements,  II,  Worldwide 
distribution,  J .  Atmos .  Sci . ,  32 : 1457-1462  . 

23.  Hofmann,  D.  J.,  and  Rosen,  J.  M.,  (1982)  Balloon-borne 

observations  of  stratospheric  aerosol  and  condensation 
nuclei  during  the  year  following  the  Mt.  St.  Helens 
Eruption,  J.  Geophys .  Res.,  87  :  11039-11061  . 

24.  Hofmann,  D.  J.,  and  Rosen,  J.  M.,  (1983)  Stratospheric 

sulfuric  acid  fraction  and  mass  estimate  for  the  1982 
volcanic  eruption  at  El  Chichon,  Geophys.  Res.  Let, 
10:313-316. 

25.  Hofmann,  D.  J.,  and  Rosen,  J.  M.,  (1984)  On  the  temporal 

variation  of  stratospheric  aerosol  size  and  mass 
during  the  first  18  months  following  the  1982  eruption 
of  El  Chichon,  J .  Geophys .  Res , ,  89 : 4883-4890 . 
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aerosols  as  having  an  average  mixture  of  75  %  and  25  % 

h2o. 

The  balloon  measurements  of  background  conditions 
yielded  average  peak  aerosol  concentrations  that  would  cor¬ 
respond  to  back  scattering  ratios  of  about  1.09  -  1.17  at 
2  6 

0.6943  ym  .  These  results  agree  well  with  the  lidar  re¬ 
sults  summarized  in  Table  1. 

Hofmann  and  Rosen  conducted  thirty  six  balloon  flights 

from  Laramie,  Wyoming  (41°  N)  to  study  the  Mt.  St.  Helens 

23 

eruption  cloud  during  the  year  following  the  eruption 
The  Mt.  St.  Helens  eruption  cloud  was  observed  as  four  dif¬ 
ferent  layers  at  different  altitudes  over  the  12  -  24  km  re¬ 
gion.  The  cloud  was  first  observed  at  about  12  -  15  km  in 
the  vicinity  of  the  jet  stream.  After  about  a  week,  a  layer 
appeared  in  the  15  -  18  km  range  over  Laramie.  This  layer 

became  the  main  layer  from  the  eruption.  Material  was  also 
injected  into  the  18  -  24  km  range  where  the  summer  strato¬ 
spheric  winds  shift  from  westerlies  to  easterlies.  The 
third  layer  was  at  18  -  20  km  where  the  winds  were  shifting 
and  the  fourth  at  20  -  24  km  where  the  winds  were  blowing 

from  the  east. 

Hofmann  and  Rosen  also  studied  the  El  Chichon  cloud  for 

n  j  o  c 

eighteen  months  following  its  erupt ion^* ' .  Those  flights 
were  from  Laramie  and  from  locations  in  southern  Texas 
(27°  N  -  29  °N) . 

26.  Pinnick,  R.  G.,  Rosen,  J.  M.f  and  Hofmann,  D.  J.,  (1976) 
Stratospheric  aerosol  measurements,  III,  Optical  model 
measurements,  J .  Atmos .  Sci . ,  33:304-314. 
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The  flights  discovered  two  major  layers  of  aerosols 
separated  by  a  very  clean  region.  The  first  was  around  the 
tropopause  and  extended  to  about  21  km.  The  second  layer 
was  centered  around  25  km  and  was  generally  about  5  km 
thick.  The  two  layers  may  have  resulted  from  two  separate 
eruptions.  The  lower  layer  may  have  resulted  from  the  ini¬ 
tial  28  March  1982  eruption  while  the  higher  one  may  have 
come  from  the  more  violent  4  April  1982  eruption. 

The  upper  layer  was  dominated  by  larger  particles  con¬ 
sisting  of  a  mixture  of  about  80  %  and  20  %  I^O  while 
the  lower  layer  consisted  of  an  aerosol  with  a  concentration 
of  60  -  65  %  sulfuric  acid.  The  difference  in  acid  concen¬ 
tration  was  due  to  the  warmer  temperatures  and  lower  water 
vapor  abundances  in  the  upper  layer. 

Table  2  shows  the  peak  aerosol  concentrations  measured 
above  20  km  from  the  two  locations  as  a  function  of  time 
following  the  eruption.  The  data  displayed  only  cover  three 
of  the  size  ranges  measured.  Data  from  a  flight  on  February 
5,  1982  are  given  as  representative  of  pre-eruption  values. 
The  differences  between  the  two  locations  result  from  the 
Texas  flights  penetrating  denser  regions  of  the  eruption 
cloud. 

The  size  distribution  determined  from  balloon  measure¬ 
ments  taken  45  days  after  the  eruption  was  bimodal  with  mode 
radii  of  0.02  and  0.7  ym.  After  about  June  1982,  the  produc¬ 
tion  of  condensation  nuclei  (e.g.,  the  curve  of  r  >  0.01  ym) 
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had  essentially  ceased  and  fallen  back  to  pre-eruption  val¬ 
ues  . 


Table  2.  Peak  Aerosol  Concentrations  Above  20  km  for  Three 
Size  Ranges  As  a  Function  of  Time  After  the  Eruption  of  El 
*  Chichon  as  Measured  from  Laramie,  Wyoming  and  Southern 

Texas^  .  The  data  from  2/5/82  are  given  as  representative 
of  pre-eruption  values 


- — - : 

PEAK  AEROSOL  CONCENTRATIONS 

DAYS  AFTER 
ERUPTION 

r  >  0 

(t 

.01,ym 
cm  J ) 

r  >  0 

(# 

.15  ym 
cm  J) 

r  > 

(« 

0. 25  ym 
cm 

Wyo 

Texas 

Wyo 

Texas 

Wyo 

Texas 

(2/5/82) 

8 

0.5 

0.1 

2 

25 

0.8 

0.1 

10 

50 

0.8 

0.1 

20 

48 

0.7 

0.1 

40 

45 

0.65 

0.13 

45 

700 

22 

16 

60 

38 

400 

0.6 

20 

0.16 

15 

73 

400 

4.0 

1.1 

80 

200 

160 

4.2 

16 

2.0 

13 

95 

45 

8.0 

6.0 

100 

38 

65 

5.0 

15 

3.0 

11 

110 

23 

2.0 

1.2 

120 

24 

30 

4.0 

12 

2.8 

10 

125 

24 

5.2 

4.0 

140 

15 

12 

4 . 0 

11 

2.8 

0.9 

152 

10 

3.5 

2.0 

160 

13 

13 

5.0 

11.5 

4.0 

0.85 

178 

10 

1.8 

0.8 

180 

11 

13.5 

2.3 

11 

1.0 

0.8 

190 

13 

10.0 

2.0 

200 

12 

14 

7.0 

10 

2.5 

0.7 

215 

11 

4.0 

3.0 

220 

10 

5.0 

3.2 

240 

10 

7.1 

5.2 

2.3  Aircraft  Measurements 

Extensive  aircraft  measurements  of  stratospheric 
aerosols  were  made  of  the  Mt.  St.  Helens  and  El  Chichon 
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eruptions.  Flight  paths  were  often  chosen  to  overlap  other 
measurements  such  as  the  balloon  measurements  of  Hofmann  and 
Rosen. 

I 

2.3.1  Size  Distribution  Measurements 

Oberbeck  et  al.  collected  particles  from  the  Mt.  St.  1 

Helens  eruption  cloud  for  a  year  using  wire  impactors  flown 
upon  the  NASA  U-2.  The  aircraft  sampled  the  stratospheric 
aerosols  at  an  altitude  of  18.3  km,  the  height  where  the 
main  aerosol  layer  eventually  formed. 

Table  3  shows  the  sampled  aerosol  particle  concentra¬ 
tions  for  two  size  ranges,  r  <  0.15  ym  and  r  >  0.15  ym,  for 
the  year  following  the  eruption.  In  the  six  months  after 
the  eruption,  the  total  volume  of  aerosols  with  radii 
greater  than  0.03  ym  increased.  The  concentration  of  parti¬ 
cles  with  radii  less  than  0.15  ym  decreased  while  the  con- 

* 

centration  of  particles  with  radii  greater  than  0.15  ym  in¬ 
creased.  This  is  consistent  with  the  dispersal  of  an  ini¬ 
tially  high  concentration  of  small  particles  and  the  growth 
of  aerosols  by  condensation.  The  aerosol  levels  returned  to 
normal  in  about  a  year. 

The  El  Chichon  eruption  was  well  studied  by  aircraft. 

One  series  of  measurements  involved  measuring  stratospheric 
size  distributions  with  "Knollenberg"  counters  .  The  mea- 

27.  Oberbeck,  V.  R.,  Farlow.  N.  H.,  Fong,  W.,  Snetsinger,  K. 

G.,  Ferry,  G.  V.  Ferry  and  Hayes,  D.  M.,  (1982)  Mount 
St.  Helens  aerosol  evolution,  Geophys.  Res,  Let., 

9:1089-1092. 

28.  Knollenberg,  R.  G.,  and  Huffmann,  D.,  (1983)  Measure¬ 

ments  of  the  aerosol  size  distributions  in  the  El 

Chichon  cloud,  Geophys.  Res.  Let.,  1£:  1025-1028 .  *» 
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surements  made  prior  to  the  eruption  revealed  relatively 
featureless  size  distributions.  The  data  over  the  range  0.1 


to  1.0  ym  revealed  particles  that  were  largely  H2S04.  Their 
number  density  varied  with  altitude  from  about  20  cm-'*  at  17 

_3 

km  to  about  1  cm  at  20  km. 


Table  3.  Aerosol  Concentrations  From  the  Mt.  St.  Helens 
Eruption  Cloud  at  18.3  km  in  two  Size  Ranges  as  Measured  by 
Wire  Impactor 


DATE 

MONTHS  AFTER 
ERUPTION 

r  >  0,15  ym 
(cm-3) 

r  <  0,15  ym 
(cm-3) 

5/20/80 

0 

2.0 

11.0 

6/25/80 

1 

1.8 

22.5 

7/16/80 

2 

2.5 

18.0 

10/29/80 

5 

3.5 

14.0 

12/2/80 

6.5 

4.0 

15.0 

12/17/80 

7 

2.5 

6.0 

4/10/81 

12 

1.5 

3.0 

Flights  following  the  eruption  suggested  two  dominant 
size  modes.  The  first  mode  consisted  of  small  particles 
(0.1  -  0.8  ym)  that  were  primarily  H2S04  and  large  particles 
(0.8  -  30  ym)  that  were  primarily  volcanic  ash.  Later 
flights  indicated  that  particles  up  to  2.0  ym  were  largely 
H2S04  although  some  of  the  particles  may  have  been  coated 
with  sulfuric  acid. 
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Oberbeck  et  al.^  also  studied  the  El  Chichon  eruption 
cloud  with  the  Ames  Wire  Impactor.  Table  4  summarizes  the 
total  aerosol  concentrations  as  measured  at  different  alti¬ 
tudes.  The  size  distribution  measured  prior  to  the  eruption 
could  be  fitted  with  a  log  normal  distribution  with  a  total 
concentration  of  between  3.4  and  4 . 3  cm  ,  a  mode  radius  of 
0.08  ym  and  a  standard  deviation  of  about  1.68.  After  the 
eruption,  the  distribution  was  best  represented  with  two  log 
normal  size  distributions,  an  enhanced  background  of  sulfu¬ 
ric  acid  particles  with  a  larger  mode  radius  and  a  sedimen¬ 
tation  mode  consisting  of  large  silicates. 


Table  4.  Total  Aerosol  Concentrations  From  the  El  Chichon 
Eruption  Clout}  at  Various  Altitudes  as  Measured  by  the  Ames 
Wire  Impactor 


DATB 

MONTHS  AFTER 
ERUPTION 

ALTITUDE 

(km) 

N_3 
(cm  3) 

10/16/81 

- 

18.3 

6.3 

5/5/82 

1 

20.7 

3.9 

7/23/82 

3.5 

20.7 

7.8 

9/23/82 

5.5 

19.8 

13.68 

11/4/82 

7 

20.7 

6.2 

11/5/82 

7 

20.7 

5.4 

12/13/82 

8 

20.7 

7.9 

29.  Oberbeck,  V.  R.,  Danielsen,  E.  F.,  Snetsinger,  K.  G.f 

Ferry,  G.  V.  Fong,  W. ,  Hayes,  D.  M. ,  (1983)  Effect  of 
the  eruption  of  El  Chichon  on  stratospheric  aerosol 
size  and  composition,  Geophys.  Res.  Let.,  10 ; 1021- 
1024. 
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30 

Wilson  et  al.  summarized  many  of  the  NASA  U-2  mea¬ 
surements  made  at  20  km  of  sub-2.5  ym  diameter  aerosols  fol¬ 
lowing  the  El  Chichon  eruption.  Their  results  indicated 
that  the  background  aerosols  were  dominated  by  particles 
with  radii  smaller  than  0.05  or  0.03  ym.  Measurements  in 
April  and  May  1982  indicated  significant  increases  in  the 
sub  0.1  ym  aerosol  concentration  that  can  only  be  explained 
by  gas-phase  reactions  leading  to  the  formation  of  secondary 
(non-eruption)  aerosols.  Measurements  in  the  same  size 
range  made  in  late  1982  showed  considerable  depletion  sug¬ 
gesting  coagulation  of  small  particles  to  form  larger  ones. 

2.3.2  Transmission  Measurements 

31 

Witteborn  et  al.  obtained  infrared  transmission  spec¬ 
tra  for  the  atmosphere  above  11  km  for  latitudes  between 
2°  S  and  50°  N.  The  measurements  were  made  in  December  1982 
when  the  stratospheric  aerosols  were  still  considerably  en¬ 
hanced  by  the  El  Chichon  eruption.  The  absorption  was  ob¬ 
tained  by  ratioing  the  transmission  at  8.5  ym,  where  the 
H2so4  absorption  is  strongest,  to  that  at  12  ym,  where  the 
atmospheric  transmission  is  nearly  unity.  The  average  ab¬ 
sorption  per  unit  air  mass  was  evaluated  at  0.019,  of  which 
0.010  was  attributable  to  the  presence  of  H2S04  aerosols. 

30.  Wilson,  J.  C. ,  Blackshear,  E.  D.,  and  Hyun,  J.  H., 

(1983)  Changes  in  the  sub-2.5  micron  diameter  aerosol 
observed  at  20  km  altitude  after  the  eruption  of  El 
Chichon,  Geophys ■  Res.  Let.,  10 : 1029-1032 . 

31.  Witteborn,  F.  C.,  O'Brien,  K.,  Crean,  H.  W. ,  Pollack,  J. 

B.,  and  Bilski,  K.  H.,  (1983)  Spectroscopic 
measurements  of  the  8-  to  13-micrometer  transmission 
of  the  upper  atmosphere  following  the  El  Chichon 
eruptions,  Geophys.  Res.  Let.,  10 ; 1009-1012 . 
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They  also  reported  that  the  ratio  of  optical  depth  at  8.5  ym 

.  32 

to  that  at  0.5  ym,  obtained  by  Dutton  and  DeLuisi  was 
about  0.14  at  20°  N. 

2.3.3  Composition  Measurements 

33 

Woods  and  Chuan  sampled  the  El  Chichon  volcanic  cloud 
using  a  quartz  crystal  microbalance  ( QCM)  cascade  impactor 
flown  on  the  NASA  U-2.  The  QCM  collected  and  classified 
particles  into  10  size  intervals  from  less  than  0.05  ym  to 
greated  than  25  ym  diameter.  The  collected  aerosols  were 
then  analyzed  for  composition. 

The  compositional  analysis  indicated  that  the  submicron 
particles  were  largely  sulfuric  acid  droplets.  This  is 
contrasted  with  Mt.  Agung  in  which  the  primary  component  was 
volcanic  ash.  Their  measurements  of  the  sulfuric  acid 
aerosols  indicated  that  they  occurred  in  a  relatively  narrow 
size  range  with  diameters  from  about  0.08  to  0.45  ym. 
Particles  above  this  range  were  primarily  lithic  and 
magmatic  materials  from  the  volcano.  Halite  particles  were 
also  found  and  were  believed  to  be  from  a  salt  dome  that  was 
located  beneath  El  Chichon. 

Patterson  and  coworkers  have  measured  the  optical  prop¬ 
erties  of  volcanic  ash  in  visble  wavelength  for  ash  from  Mt. 

32.  Dutton,  E.,  and  DeLuisi,  J.,  (1983)  Extinction  of  solar 

radiation  by  the  El  Chichon  dust  cloud,  Geophys .  Res. 
Let . ,  10: 1013-1016. 

33.  Woods,  D.  C.,  and  Chuan,  R.  L.,  (1983)  Size-specific 

compositions  of  aerosols  in  the  El  Chichon  volcanic 
cloud,  Geophys.  Res.  Let.,  10:1041-1044. 
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St.  Helens  and  El  Chichon  .  Their  studies  indicate  that 
the  imaginary  component  of  the  refractive  index  of  the  ash 
varies  from  volcano  to  volcano.  The  Mt.  St.  Helens  ash  had 
an  imaginary  component  that  decreased  from  about  0.01  at 
0.3  pm  to  about  0.0015  at  0.7  pm.  The  corresponding  values 
for  the  El  Chichon  ash  were  nearly  constant  at  0.001  over 
the  same  wavelength  range. 

2.4  Si2e  Distribution  Studies  Using  Satellite  Data 

A  new  series  of  satellites  have  been  launched  that  of¬ 
fer  new  details  of  the  background  and  perturbed  strato¬ 
sphere.  These  satellites  have  include  the  Stratospheric 
Aerosol  and  Gas  Experiment  (SAGE  I  &  II)  and  the  Strato¬ 
spheric  Aerosol  Measurement  (SAM  I  &  II)  sensors.  Although 
not  explicitly  designed  to  study  stratospheric  aerosols,  the 
Solar  Mesosphere  Explorer  spacecraft  was  also  able  to  study 
the  El  Chichon  eruption. 

The  SAM  II  sensor  has  a  single  pass  band  centered  at 
1.0  pm.  Gaseous  absorption  is  negligible  at  this  wave¬ 
length,  so  the  measured  extinction  is  due  to  molecular  and 
aerosol  scattering.  The  molecular  scattering  is  filtered 

out  during  data  processing  leaving  only  the  aerosol 

3  6 

component.  Russell  et  al.  performed  a  comparison  of  SAM 

34.  Patterson,  E.  M.  (1981)  Measurements  of  the  imaginary 

part  of  the  refractive  index  between  300  and  700  nano¬ 
meters  for  Mount  St.  Helens  ash,  Science,  211 ; 836-838 . 

35.  Patterson,  E.  M.,  Pollard,  C.  0.,  and  Galindo,  I., 

(1983)  Optical  properties  of  the  ash  from  El  Chichon, 
Geophys.  Res.  Let. ,  10:317-320. 

36.  Russell,  P.  B.,  McCormick,  M.  P.,Swissler,  T.  J.  Chu, 

W.P.,  Livingston,  J.  M.,  Fuller,  W.  H  ,  Rosen,  J. 

M.,  Hofmann,  D.  J.,  McMaster,  L.  R.,  Woods,  D.  C.  and 
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II  extinction  profiles  against  profiles  inferred  from  bal¬ 
loon-borne  dustsondes  and  an  airborne  lidar.  Their  results 
indicated  that  the  different  sensors  agreed  within  their  re¬ 
spective  measurement  and  conversion  uncerta int ies . 

The  SAGE  I  sensor  system  consists  of  four  wavelengths 
with  two  of  them,  centered  at  0.45  and  1.0  um,  designed  for 
aerosol  studies.  The  two  channels,  used  together,  provide 

extinction  profiles  and  size  distribution  data. 

37 

Lenoble  et  al.  used  93  SAGE  I  profiles  of  the  ex¬ 
tinction  coefficients  at  1.0  and  0.45  ym  from  July  1980  to 
study  the  variation  of  the  aerosol  size  distributions  in  the 
Mount  St.  Helens  aerosol  layer.  The  measurements  were  made 
between  50  and  70°  N.  Further  measurements  from  May  to 
November  1980  showed  a  low  layer  cloud  of  primarily  large 
particles  and  a  higher  altitude  cloud  of  small  particles. 
The  size  distributions  within  the  clouds  were  assumed  to  be 
log  normal  and  the  clouds  were  assumed  to  be  made  up  of  75  % 
H2S04  droplets  in  solution  with  H20.  Using  Mie  theory,  the 
mode  radii  of  the  layers  were  inverted.  The  low  layer  cloud 
was  found  to  have  an  average  logarithmic  mode  radius  of 
0.20  urn  while  the  high  altitude  cloud  had  a  logarithmic  mode 
radius  of  0.06  ym.  Logarithmic  mode  radii  corresponding  to 

Pepin,  T.  ( 1981 )  Satellite  and  correlative 
measurements  of  the  stratospheric  aerosol.  II: 
Comparison  of  measurements  made  by  SAM  II,  dust¬ 
sondes  and  an  airborne  Lidar,  J.  Atmos.  Sci.,  38: 
1295-1312. 

37.  Lenoble.  J.,  Pruvost,  P.  and  Brogniez,  C.  (1984)  SAGE 

satellite  observations  of  stratospheric  aerosols  from 
Mount  St.  Helens  eruption:  A  two-wavelength  analysis, 

J.  Geophys.  Res.,  89 : 11666-11676 . 
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individual  orbits  were  given  for  approximately  0.2  degree 
intervals  from  50°  to  65°  N  latitudes. 

Lati tudinally  averaged  results  were  given  as  a  function 
of  time  that  demonstrate  the  temporal  variation  of  the  size 
distributions  as  the  larger  particles  settle  out.  Those  re¬ 
sults  are  given  in  Table  5. 


Table  5.  Lat i tudinal ly  Averaged  Mode  Radii,  in  ym,  for  a 
Log  Normal  Size  Distribution  as  a  Function  of  Time  Between 
55°  N  and  65°  N  Latitudes  Following  the  Mount  St.  Helens 
Eruption-*  .  These  values  assume  a  standard  deviation  of 
1.604 


June  23  to  July  20  to  Aug  28  to  Sept  28  to 
July  20  Aug  28  Sept  28  Oct  31 


iininn 

(km) 

55°  65° 

55° 

65° 

55° 

65° 

55° 

65° 

18-20 

22-23 

---  0.151 

-  0.084 

0.142 

0.090 

0.146 

0.073 

0.142 

0.100 

0.137 

0.100 

0.142 

0.106 

0.142 

0.103 

O  Q 

Brogniez  and  Lenoble  have  also  used  SAGE  I  profiles 
to  study  background  stratospheric  aerosol  conditions.  From 
the  tropopause  to  about  10  km  above  the  tropopause,  the 
background  extinction  coefficient  at  1.0  p  is  nearly  con¬ 
stant  at  all  latitudes  and  seasons  at  1.20  x  10-<*  km.  Above 
this  altitude,  the  background  extinction  coefficient  profile 
can  be  fitted  with  an  exponential  with  an  average  scale 
height  of  3.2  km.  The  altitude  where  the  profiles  change 

38.  Brogniez,  C.,  and  Lenoble,  J.,  (1987)  Modeling  of  the 
stratospheric  background  aerosols  from  zonally 
averaged  SAGE  profiles,  J.  Geophys.  Res.,  92 : 3051- 
3060. 
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from  constant  to  an  exponential  decrease  is  a  function  of 
latitude  and  season.  The  transition  altitude  is  highest  for 
high  latitude  winters  and  lowest  for  high  latitude  summers. 

Their  values  can  be  compared  to  those  in  the  (SRA) 
model  .  The  SRA  model  gives  an  extinction  coefficient  of 
0.68  x  10-4  km-1  for  1.0  ym  for  the  altitudes  12  -  20  km. 
Above  20  km,  the  decrease  in  extinction  coefficient  is 
linear . 
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3 .  METHODOLOGY 


3.1  Calculating  the  Temperature  Dependence  of  the  Index  of 
Refraction  of  ^SO^  Droplets 

The  temperature  dependence  of  sulfuric  acid  droplets 
for  wavelengths  far  from  absorption  features  can  be  calcu¬ 
lated  using  the  Lorentz-Lorenz  relationship.  The  Lorentz- 
Lorenz  relationship  relates  the  complex  index  of  refraction, 
m,  to  the  molecular  polarizability  per  unit  volume 


P  = 


3V 


4  TT 


(m  -  1) 
(m2  +  2) 


(2) 


where  p  is  the  density.  The  Lorentz-Lorenz  equation  is  a 

form  of  the  Clausius-Mossotti  equation  that  established  that 

( e-  l)/(e+  2)  should  be  proportional  to  the  density  of  a 

substance,  where  e  is  the  dielectric  constant.  Eq.  (2)  is 

2  39 

obtained  by  setting  e  =  m  .  According  to  Jackson  and 

40 

Mayer  and  Mayer  ,  the  proportionality  to  density  holds  for 

gases  but  is  only  an  approximation  for  liquids.  Pinkley  and 
g 

Williams  have  observed  experimentally  that  the  Lorentz- 
Lorenz  relationship  appears  to  be  valid  away  from  absorption 
bands,  but  fails  near  strong  bands  for  liquids. 

Under  the  assumption  that  the  relationship  is  valid  for 
a  solution  of  I^SC^  and  water,  the  Lorentz-Lorenz  relation¬ 
ship  can  be  expressed  as 


39.  Jackson,  J.  D.  (1975)  Classical  Electrodynamics,  John 

Wiley  &  Sons,  2  nd  Edition,  New  York,  p  155. 

40.  Mayer,  J,  E.,  and  Mayer,  M.  G.  (1940)  Statistical 

Mechanics,  John  Wiley  &  Sons,  New  York,  p  340. 
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1) 


=  A  p iT) 


(3) 


(m2  - 


(m  +  2) 


where  A  is  a  constant  and p  (T)  is  the  temperature  dependent 
solution  density.  Letting  be  the  index  of  refrac¬ 
tion  at  a  temperature  and  m2  be  the  index  of  refraction 

at  temperature  T2,  then 


-  1) 
+  2) 


=  £-1 
p  1 


(m2  -  1) 
(m2  +  2) 


(4) 


or,  with  a  little  algebra 


m2  = 


'  1  +  2D  "1 

.  1  -  D  J 


1/2 


(5) 


where  D  is  given  by 

D  =  £—2. 
P  1 


(m2  -  1) 
(m2  +  2) 


(6) 


The  temperature  dependent  density  data  of  a  75  %  solution  of 
41 

H2SC>4  and  water  used  in  the  calculations  are  shown  in  Fig¬ 
ure  1 . 

For  wavelength  regions  in  the  vicinity  of  absorption 

features,  the  Lorentz-Lorenz  relationship  is  not  valid.  In 

those  regions,  a  linear  interpolation  and  extrapolation  of 

the  available  data  was  performed.  For  the  imaginary 

components,  the  interpolation  was  performed  in  terms  of  the 

logarithm  of  the  value.  For  wavelengths  above  3.0  ym,  the 

41.  Roth,  W.  A.  and  Scheel,K.,  eds.  (1923)  Landolt-Born- 
stein  Physikalisch-Chemische  Tabellen,  Verlag  von 
Julius,  Springer,  Berlin,  397-398. 


22 


75%  H-SO.  and 
The  flashed 

indicates  the  extrapolated  values 


interpolation 

procedure 

was  performed 

using  the 

data 

discussed  below  at 

300  K 

and  the  data 

of 

Pinkley 

and 

9 

Williams  at 

250  K. 

The 

earlier  values 

of 

Shettle 

and 

Fenn1,  based 

_  7 

on  an 

average  of  the  data 

of 

Palmer 

and 

7  c  & 

Williams  and  Remsberg  '  were  used  for  wavelengths  up  to 


A  0 

25  nm.  Between  20  and  50  ym,  the  data  of  Jones  were  used. 

Between  50  and  300  ym,  an  interpolation  of  qualitative 

43 

results  from  the  work  of  Rossow  and  Sagan  were  used  to 


42.  Jones,  A.  D.  (1976)  Optical  constants  of  sulfuric  acid 

in  the  far  infrared,  J.  Quant.  Spectrosc.  Radiat. 
Transfer,  16:1017-1019. 

43.  Rossow,  W.  B.,  and  Sagan,  C.,  (1975)  Microwave  boundary 

conditions  on  the  atmosphere  and  clouds  of  Venus,  J. 
Atmos .  Sci . ,  32:1164-1176. 


25%  H2 
portio 


Figure  1.  The  Density  of  a  Mixture. of 
as  a  Function  of  Temperature41. 
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extend  the  values  for  both  temperatures.  Those  results 
assumed  that  the  imaginary  part  of  the  index  of  refraction 
reached  a  maximum  at  150  cm  1  and  decreased  linearly  with 
frequency  to  about  10  cm  1 .  The  real  part  was  assumed  to 
varying  linearly  with  frequency  between  40  and  1000  ym, 
reaching  about  2.1  at  1000  ym. 

3.2  Proposed  Size  Distribution 

The  log  normal  size  distribution  parameters  suggested 
by  the  World  Meteorological  Organization's  Experts  Meeting10 
will  be  used.  These  parameters  are  consistent  with  those 
based  on  the  review  of  various  measurements  discussed  in 
Chapter  2.  The  distribution  is  given  by 


The  extinction  coefficient  values  produced  by  (7)  have  been 

-4  -1 

scaled  to  give  1.2  x  10  km  at  1 . 0  ym,  the  average  SAGE  I 

38 

value  reported  by  Brogniez  and  Lenoble  .  The  resultant 

_  3 

values  are  a  total  number  density,  Nq,  of  9  cm  ,  a  mode 
radius,  R,  of  0.06951  ym  and  a  standard  deviation,  0,  of 
1.86.  Figure  2  compares  the  modified  gamma  function  used  by 
Shettle  and  Fenn1  and  the  proposed  log  normal  size 
distribution . 
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Figure  2.  Comparison  of  the  Modified  Gamma  Function  Size 
Distribution  Used  in  Shettle  and  Fenn1  and  the  Proposed  Log 
Normal  Size  Distribution 
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4 .  RESULTS 


4.1  Indices  of  Refraction 

Calculations  for  temperature  dependent  indices  of  re¬ 
fraction  were  performed  for  temperatures  from  150  to  290  K. 
The  calculations  were  done  for  68  wavelengths  between  0.2 
and  300  pm.  The  Lorentz-Lorenz  relationship  was  used  for 
wavelengths  below  3.0  vm  with  data  obtained  at  300  K  used  as 
reference  .  For  wavelengths  above  3.0  pm,  values  were  in¬ 
terpolated  as  described  in  Section  3.1.  The  values  are 
listed  in  Table  6. 

Figure  3  shows  the  temperature  dependence  of  the  real 

component  of  m  at  0.55  pm.  At  this  wavelength,  the  imagi- 

_  8 

nary  component  is  on  the  order  of  10  and  is  nearly  temper¬ 
ature  independent.  The  inflection  points  seen  near  160,  220 
and  280  K  are  a  result  of  truncating  the  refractive  index 
values  to  three  significant  figures.  Figures  4  (a)  and  (b) 
show,  respectively,  the  real  and  imaginary  components  at 
10  pm  as  a  function  of  temperature. 

Indices  of  refraction  at  215  K  will  be  taken  as  repre¬ 
sentative  of  stratospheric  temperatures  and  used  to  deter¬ 
mine  the  extinction  coefficients  for  the  proposed  background 
stratospheric  aerosol  model.  The  indices  of  refraction  at 
215  K  are  listed  in  Table  6  and  the  real  and  imaginary  parts 
shown,  respectively,  in  Figures  5  (a)  and  (b)  as  a  function 
of  wavelength.  The  values  at  300  K,  at  250  K,  and  at  150K 
are  also  shown  for  comparison. 


Figure  3.  The  Real  Part  of  the  Index  of  Refraction  at 
0.55  ym  of  a  Mixture  of  75%  H2S04  and  25%  H20  as  a  Function 
of  Temperature.  The  inflection  points  that  appear  near  160, 
220,  and  280  K  are  a  result  of  truncating  the  refractive 
index  values  to  three  significant  figures 
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IMAGINARY  PART 


150  200  250  300 

TEMPERATURE  (K) 


Figure  4.  The  (a)  Real  and  (b)  Imaginary  Parts  of  the  Index 
of  Refraction  at  10  ym  of  a  Mixture  of  75%  H-SO.  as  a 
Function  of  Temperature 
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Figure  5.  The  (a)  Real  and  (b)  Imaginary  Parts  of  the  Index 
of  Refraction  for  a  75%  Sulfuric  Acid  Solution  at  300,  250, 
215  K  and  150  K  as  a  Function  of  Wavelength 


Table  6.  Real,  n  ,  and  Imaginary,  n.,  Indices  of  Refraction 
for  a  75  %  Solution  of  Sulfuric  Acici  Droplets  at  215  K.  The 
values  in  parentheses  indicate  the  power  of  10  by  which  n. 
should  be  multiplied  1 


f - - -  - -  - -  - -  •  -  1 

LAMBDA 

(ym) 

nr 

ni 

LAMBDA 

(urn) 

nr 

ni 

0.2000 

1 .526 

-1 . 07 ( -8 ) 

8.5000 

1 .457 

-0 . 840 

0.2500 

1.512 

-1 . 07 ( -8  ) 

8.7000 

1.753 

-0.817 

0.3000 

1.496 

-1 . 07 ( -8 ) 

9.0000 

1.767 

-0.604 

0 .3371 

1.484 

-1 . 07 ( -8  ) 

9 .2000 

1.704 

-0.557 

0.4000 

1 .464 

-1 . 07 ( - 8  ) 

9 .5000 

1.791 

-0.736 

0 . 4880 

1 .456 

-1 . 0  7 ( -  8  ) 

9 .8000 

2.128 

-0.465 

0 . 5145 

1.454 

-1 . 07 ( -8 ) 

10.0000 

2.094 

-0.306 

0 . 5500 

1 .454 

-1 . 07 ( -8 ) 

10 . 5910 

1.805 

-0.210 

0 .6328 

1.452 

-1 . 56 ( -8 ) 

11.0000 

1.823 

-0.433 

0 .6943 

1.452 

-2 . 12 ( -8 ) 

11.5000 

2.024 

-0.199 

0.8600 

1.448 

-1 . 90 ( -7 ) 

12 . 5000 

1.808 

-0.110 

1.0600 

1.443 

-1 . 60 ( -6 ) 

13.0000 

1.790 

-0.112 

1.3000 

1.432 

-1 . 06 ( -5  ) 

14.0000 

1.739 

-0.134 

1 .5360 

1.425 

-1 .46(-4  ) 

14 . 8000 

1.663 

-0.163 

1.8000 

1.411 

-  5 . 8  5  (  -  4  ) 

15.0000 

1.656 

-0.181 

2.0000 

1.405 

-0.00134 

16.4000 

1.692 

-0.463 

2 .2500 

1.390 

-0.00192 

17.2000 

2.003 

-0.616 

2.5000 

1.362 

-0.00399 

18.0000 

2.207 

-0.166 

2.7000 

1 . 319 

-0.00604 

18.5000 

2.041 

-0.0940 

3.0000 

1 .288 

-0.0875 

20.0000 

1.895 

-0.0558 

3.2000 

1.292 

-0.144 

21 . 3000 

1.809 

-0.146 

3.3923 

1.366 

-0.178 

22.5000 

1.942 

-0.187 

3.5000 

1.383 

-0.179 

25.0000 

1.993 

-0.0274 

3.7500 

1.420 

-0.165 

27.9000 

1 . 916 

-0.0221 

4.0000 

1.435 

-0.151 

30.0000 

1.866 

-0.0259 

4 . 5000 

1.427 

-0.156 

35 . 0000 

1.856 

-0.0315 

5.0000 

1.411 

-0.137 

40.0000 

2.093 

-0.131 

5.5000 

1.376 

-0.195 

50.0000 

2.163 

-0.238 

6.0000 

1.485 

-0.195 

60.0000 

2.163 

-0.263 

6 . 2000 

1.485 

-0.147 

80.0000 

2.156 

-0.276 

6.5000 

1.421 

-0.0841 

100.0000 

2.173 

-0.261 

7.2000 

1 . 235 

-0.166 

150.0000 

2.183 

-0.172 

7.9000 

1.148 

-0.461 

200.0000 

2.176 

-0.107 

8.2000 

1.220 

-0.709 

300.0000 

2.179 

-0.0386 

4.2  The  Proposed  Background  Stratospheric  Aerosol  Model 

Figure  6  shows  the  extinction,  scattering  and  absorp¬ 
tion  coefficients  of  the  proposed  background  stratospheric 
aerosol  model.  Table  7  lists  the  values  along  with  the 
asymmetry  parameters.  Figures  7  (a)  to  (d)  compare  the 
respective  parameters  against  those  developed  by  Shettle  and 
Fenn1  and  used  in  LOWTRAN  6^. 

The  extinction  coefficients  for  the  proposed  model  are 
smaller  than  those  used  of  Shettle  and  Fenn.  The  differ¬ 
ences  are  greatest  for  wavelengths  longer  than  10  ym.  The 
differences  are  almost  entirely  due  to  changes  in  the  ab¬ 
sorption  coefficient,  as  shown  in  Figure  7  (c).  The  differ¬ 
ences  in  the  scattering  coefficients  and  asymmetry  parameter 
are  slight. 

The  changes  seen  in  the  proposed  formulation  are  due 
almost  entirely  to  the  new  indices  of  refraction.  Figure  8 
compares  the  extinction  coefficients  as  a  function  of  wave¬ 
length  of  the  proposed  model  against  those  calculated  using 
the  proposed  indices  of  refraction  and  the  old  modified 
gamma  size  distribution  and  against  those  calculated  using 
the  indices  of  refraction  at  300  K  and  the  proposed  log  nor¬ 
mal  size  distribution.  As  shown,  the  change  in  indices  of 
refraction  is  responsible  for  the  majority  of  the  differ¬ 
ence  . 

4.2.1  Impact  of  Choice  of  Size  Distribution 

Based  on  our  review  of  the  literature  (see  Chapter  2) 
the  size  distribution  that  was  chosen  follows  rec- 
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Table  7.  Radiative  Parameters  for  the  Proposed  Background 
Stratospheric  Aerosol  Model  at  215  K.  The  numbers  in 
parentheses  are  powers  of  ten 


< — — - - - - 

LAMBDA 

( y«) 

INDEX  OF 
REFRACTION 

nr  ni 

EXTINCTION 
( kja  1 ) 

SCATTERING 

ABSORPTION 
Oun  l) 

SINGLE 

SCATTERING 

ALBEDO 

ASYMMETRY 

PARAMETER 

0.2000 

1.526  0.0000 

8.50421-4) 

8.50421-4) 

9.97741-11  ) 

1.000 

0.6749 

0 .2500 

1.512  0.0000 

8 . 2526 ( -4 ) 

8.25261-4) 

7.76181-11) 

1  .000 

0.6850 

0.3000 

1.496  0.0000 

7.64541-4) 

7.64541-4) 

6.12601-11  ) 

1.000 

0.6943 

0.3371 

1.48'  0.0000 

7.0676(-4) 

7.06761-4 ) 

5.23561-11) 

1.000 

0.6991 

0.4000 

1.464  0.0000 

5.99641-4) 

5.99641-4) 

4.23641-11 ) 

1.000 

0.7038 

0.4880 

1.456  0.0000 

4.79101-4) 

4.79101-4) 

3.29201-11) 

1.000 

0.6944 

0.5145 

1.454  0.0000 

4.46941-4) 

4.46941-4) 

3.06421-11) 

1.000 

0.6910 

0.5500 

1.454  0.0000 

4.09801-4 ! 

4.09801-4  ) 

2.81601-11) 

1 .000 

0.6846 

0.6328 

1.452  0.0000 

3.33181-4) 

3.33181-4) 

3.42441-11) 

1.000 

0.6694 

0.6943 

1.452  0.0000 

2.87501-4) 

2.87501-4) 

4.11961-11  ) 

1 .000 

0.6572 

0.8600 

1.448  0.0000 

1.93521-4) 

1.93521-4) 

2.75741-10) 

1.000 

0.6251 

1.0600 

1.443  0.0000 

1.23641-4) 

1  .23631-4  ) 

1.73751-9) 

1.000 

0.5861 

1.3000 

1.432  0.0000 

7.37061-5) 

7.36961-5) 

8.63781-9) 

1.000 

0.5405 

1.5360 

1.425  -0.0001 

4.68801-5) 

4.67861-5) 

9.45601-8) 

0.998 

0 .4965 

1.8000 

1.411  -0.0006 

2.87321-5) 

2.84261-5) 
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Figure  6.  Extinction,  Scattering  and  Absorption 
Coefficients  as  a  Function  of  Wavelength  for  the  Proposed 
Background  Stratospheric  Aerosol  Model 
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Figure  7.  Comparison  of  the  Proposed  Background 
Stratospheric  Aerosol  Model  and  that  in  Shettle  and  Fenn  (a) 
Extinction  Coefficients,  (b)  Scattering  Coefficients,  (c) 
Absorption  Coefficients,  and  (d)  Asymmetry  Parameter 
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Figure  8.  Extinction  Coefficients  as  a  Function  of 
Wavelength  for  the  Proposed  Model  (Solid  Line),  Calculated 
Using  the  Proposed  Indices  of  Refraction  and  the  Modified 
Gamma  Size  Distribution  (Dotted  Line)  and  Calculated  Using 
the  Indices  of  Refraction  at  300  K  and  the  Proposed  Log 
Normal  Size  Distribution 
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ommend^tions  from  the  meeting  of  experts  held  by  the  World 
Meteorological  Organization^.  Other  sets  of  log  normal  pa¬ 
rameters  have  been  suggested  by  other  researchers  and  Fig¬ 
ure  9  compares  the  extinction  coefficients  as  a  function  of 
wavelength  calculated  with  the  proposed  size  distribution 
against  those  calculated  by  other  suggested  size  distribu¬ 
tions.  The  indices  of  refraction  at  215  K  were  used  in  all 
cases . 


Figure  9.  Extinction  Coefficients  as  a  Function  of  Wave¬ 
length  for  the  Proposed  Model  and  for  Three  Alternative  Size 
Distributions 

The  values  of  the  parameters  for  the  log  normal  distri¬ 
bution  are  listed  in  Table  8.  The  particle  number  densities 

were  scaled  so  as  to  produce  an  extinction  coefficient  of 
-4  -1 

1.2  x  10  km  at  1.0  ym.  The  curves  labled  Lenoble  1  and 

3  ft 

2  are  from  the  work  of  Brogniez  and  Lenoble  and  the  curve 

2  6 

labled  Pinnick  is  from  Pinnick  et  al.  .  The  curves  are  es- 


Table  8.  Values  of  Log  Normal  Size  Distribution  Parameters 
in  Eq.  (7)  for  the  Different  Stratospheric  Models  Compared 
in  Figure  9.  The  N  values  have  beennormalized  to  produce 
an  extinction  coefficient  of  1.2  x  10  km  at  the  wave- 
•  length  1.0  ym 


1  _ _ _ _ _ 1 

CURVE 

(«  cm 

R 

(yrn) 

0 

Proposed 

9.0 

0.0695 

1.86 

Lenoble  1 

3.3 

0.12 

1.604 

Lenoble  2 

11.7 

0.06 

1.919 

Pinnick 

7.4 

0.0725 

1.86 

sentially  the  same,  as  shown  in  Figure  10  which  shows  the 
ratios  of  the  extinction  coefficients  shown  in  Figure  9  to 
the  proposed  values,  leading  one  to  draw  the  conclusion  that 
the  choice  of  parameters  for  the  log  normal  distribution  has 
little  impact  on  the  resultant  extinction  coefficients. 

4.2.2  Impact  of  Choice  of  Stratospheric  Temperature 

The  value  of  215  K  was  chosen  to  be  representative  of 
stratospheric  conditions  in  an  annually  averaged,  mid-lati¬ 
tude  atmosphere.  The  impact  of  using  these  values  for  other 
latitudes  or  seasons  can  be  demonstrated  in  Figure  11  which 
shows  extinction  coefficients  at  250  K  and  180  K.  The  val¬ 
ues  at  250  K  are  greater  for  the  far  infrared  wavelengths  as 
a  result  of  the  larger  imaginary  component  (see  Appendix  A). 
Conversely,  the  lower  imaginary  components  at  180  K  are  re¬ 
flected  as  lower  extinction  coefficients  at  the  far  infrared 
wavelengths.  The  values  of  the  attenuation  coefficients  for 
the  extreme  stratospheric  temperatures  of  180  and  250  K  are 
tabulated  in  the  Appendix  for  reference. 
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Figure  10.  Ratios  of  the  Extinction  Coefficients  Shown  in 
Figure  9  to  the  Proposed  Formulation  as  a  Function  of  Wave¬ 
length 
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Figure  11.  Extinction  Coefficients  as  a  Function  of  Wave 
length  for  the  Proposed  Model  (and  Assumed  Stratospheric 
Temperature  of  215  K)  and  With  Indices  of  Refraction 
Calculated  for  180  K  and  250  K 
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4.3  Enhanced  Aerosols  After  a  Volcanic  Eruption 

A  volcanic  eruption  will  increase  the  numbers  and  size 
distribution  of  sulfuric  acid  particles  as  a  result  of  S02 
injections  into  the  stratosphere.  The  specific  impact  on 
aerosol  loadings  in  terms  of  size  distribution  and  mixture 
of  H2S04  and  water  will  vary  from  eruption  to  eruption. 

Figure  12  shows  the  extinction  coefficients  for  strato¬ 
spheric  aerosols  that  have  been  "enhanced"  by  the  chemical 
production  and  sedimentation  of  sulfuric  acid  particles. 

The  calculations  were  based  on  the  sum  of  two  log  normal 

2  9 

size  distributions  based  upon  the  work  of  Oberbeck  et  al. 
following  the  eruption  of  El  Chichon.  The  first  distribu¬ 
tion  represents  an  enhanced  background  and  has  a  particle 


WAVELENGTH  (/zm) 


Figure  12.  Extinction  Coefficients  as  a  Function  of  Wave¬ 
length  for  a  Postulated  Enhanced  Stratospheric  Aerosol  Layer 


41 


density  of  6  cm’^,  mode  radius  of  0.14  ym  and  standard 
deviation  of  1.72  (e.g.  Eq.  (7)).  The  second  distribution 
results  from  the  sedimentation  of  smaller  particles  and  is 

_  3 

represented  by  a  particle  density  of  1 . 5  cm  ,  a  mode  radius 
of  0.54  ym  and  a  standard  deviation  of  1.22.  The  calcula¬ 
tions  were  done  with  the  proposed  indices  of  refraction  at 
215  K  and  those  at  300  K. 

Both  curves  show  a  nearly  flat  response  as  a  function 
of  wavelength  up  to  about  1.0  ym.  Significant  differences 
between  the  curves  do  not  begin  to  appear  until  beyond  20  ym 
as  a  result  of  the  differences  in  the  imaginary  component  of 
the  indices  of  refraction. 


4 
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5.  CONCLUSIONS  AND  RECOMMENDATIONS 
5.1  Conclusions 

0  A  new  background  stratospheric  aerosol  model  for  use  in 

LOWTRAN  7^  has  been  developed.  The  new  model  is  intended 

►  to  replace  the  background  stratospheric  aerosol  model 

2 

formulation  currently  in  use  in  LOWTRAN  6  that  is  based  on 
indices  of  refraction  for  sulfuric  acid  droplets  at  300  K. 
The  new  model  includes  indices  of  refraction  for  a  mixture 
of  75  %  and  25  %  at  temperatures  representative  of 

a  mean  stratosphere  (215  K)  and  a  log  normal  size  distribu¬ 
tion  that  is  consistent  with  observations.  Indices  of  re¬ 
fraction  at  temperatures  between  150  and  300  K  have  also 
been  calculated. 

*  The  extinction  coefficients  for  the  proposed  model  are 
smaller  than  those  used  in  LOWTRAN  6.  The  differences  are 
greatest  for  wavelengths  longer  than  10  ym.  The  differences 
are  almost  entirely  due  to  changes  in  the  absorption  coeffi¬ 
cient  as  a  result  of  smaller  imaginary  components  of  the  in¬ 
dex  of  refraction  at  these  wavelength. 

5.1.1  Impact  of  Choice  of  Size  Distribution 
There  is  little  impact  from  the  change  from  a  modified 
gamma  to  a  log  normal  size  distribution.  In  addition,  the 
choice  of  parameters  for  the  log  normal  size  distribution 
has  little  impact  on  the  extinction  coefficients. 

*  44.  Kneizys,  F.X.,  Shettle,  E.P.,  Abreu,  L.W.,  Chetwynd, 

J.H.,  Anderson,  G.P.,  Gallery,  W.O.,  Selby,  J.E.A., 
and  Clough,  S.A.  (1988)  Users  Guide  to  LOWTRAN  7, 
AFGL-TR-88-0177 ,  Air  Force  Geophysics  Laboratory, 
Hanscom  AFB ,  Massachusetts,  16  August  1988.  ADA206773 
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5.1.2  Impact  of  Choice  of  Representative  Stratospheric 
Temperature 

The  value  of  215  K  was  chosen  to  be  representative  of 
stratospheric  conditions  in  an  annually  averaged,  mid-lati¬ 
tude  atmosphere.  Examples  of  extinction  coefficients  deter¬ 
mined  for  indices  of  refraction  representative  of  extreme 
stratospheric  conditions  have  also  been  presented. 

5.2  Recommendations  for  Future  Research 

A  single  set  of  radiative  parameters  for  the  proposed 
background  stratospheric  aerosol  model  have  been  calculated 
in  order  to  simplify  its  inclusion  in  LOWTRAN  7  .  While 
the  temperature  chosen  is  more  representative  of  tempera¬ 
tures  found  in  the  stratosphere,  it  would  not  correctly  rep¬ 
resent  conditions  found  at  low  or  high  latitudes  or  be 
representative  of  seasonal  changes. 

Values  of  the  indices  of  refraction  for  at  tem¬ 
peratures  representative  of  the  full  range  of  stratospheric 
temperatures  were  also  calculated.  These  indices  of  refrac¬ 
tion  could  be  coupled  to  a  temperature  dependent  radiation 
model  in  which  the  radiative  parameters  are  coupled  to  a 
mean  temperature  inputed  for  the  stratosphere  by  the  user. 
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Appendix  A 


Temperature  Dependent  Indices  of  Refraction  and  Radiative 
Properties  for  the  Background  Stratospheric  Model 

Table  A-l  lists  the  temperature  dependent  indices  of 
refraction  for  sulfuric  acid  and  water  for  the  temperature 
range  from  150  to  300  K.  Table  A-2  lists  the  radiative 
properties  for  the  background  stratospheric  model  at  180  K 
and  Table  A-3  lists  the  properties  at  250  K. 


A-l 


Table  A-l.  Temperature  Dependent  Real,  nr,  and  Imaginary, 
n^,  Indices  of  Refraction  for  a  75%  Solution  of  H-SO^  and 
Water.  The  numbers  in  parentheses  are  powers  of  ten 


LAMBDA 
(  Vm) 

300  K 

n  n- 

r  l 

290  K 

nr  ni 

280  K 

nr  ni 

270  K 

nr  ni 

0.2000 

1.498  -1 . 00 ( -8 ) 

1.501  -1.0K-8) 

1.505  -1.02(-8) 

1.508  -1 . 02 ( -8 ) 

0.2500 

1.484  -1 . 00 ( -a ) 

1.487  -1.0K-8) 

1.490  -1 . 02 { -8  ) 

1.494  -1 . 02 ( -8 ) 

0.3000 

1.469  -1 . 00 ( -8 ) 

1.472  -1.011-8) 

1.475  -1 . 02 ( -8 ) 

1.478  -1 . 02 ( -8 ) 

0.3371 

1.459  -1 . 00 ( -8 ) 

1.462  -1.0K-8) 

1.465  -1 . 02 ( -8  ) 

1.468  -1 . 02 ( -8 ) 

0.4000 

1.440  -1 . 00 ( -8 ) 

1.443  -l.Olt-8) 

1.446  -1 . 02  ( -8 ) 

1.449  -1 . 02 ( -8 ) 

0.4880 

1.432  -1 . 00 ( -8 ) 

1.435  -1.0K-8) 

1.438  -1 . 02 ( -8 ) 

1.441  -1 . 02 ( -8 ) 

0.5145 

1.431  -1 . 00 ( -8 ) 

1.434  -1.011-8) 

1.437  -1 . 02 ( -8  ) 

1.439  -1 . 02 ( -8 ) 

0.5500 

1.430  -1 . 00 ( -8 ) 

1.433  -1.011-8) 

1.435  -1 . 02 ( -8 ) 

1.438  -1 . 02 ( -8 ) 

0.6328 

1.429  -1 . 47 ( -8 ) 

1.432  -1.481-8) 

1.434  -1 . 49 ( -8 ) 

1.437  -1 . 50 ( -8 ) 

0.6943 

1.428  -1 . 99 ( -8  ) 

1.431  -2.001-8) 

1.433  -2.021-8) 

1.436  -2 . 04 ( -8  ) 

0.8600 

1.425  -1 . 79 ( -7 ) 

1.428  -1.801-7) 

1.430  -1.821-7) 

1.433  -1 . 83 ( -7 ) 

1.0600 

1.420  -1 . 50 ( -6 ) 

1.423  -1.511-6) 

1.425  -1.521-6) 

1.428  -1 . 53 ( -6 ) 

1.3000 

1.410  -1 . 00 ( -5 ) 

1.413  -1.011-5) 

1.415  -1.011-5) 

1.418  -1 . 02 ( -5 ) 

1.5360 

1.403  -1 . 37 ( -4 ) 

1.405  -1.381-4) 

1.408  -1.391-4) 

1.411  -1 . 40 ( -4 ) 

1.8000 

1.390  -5 . 50 ( -4 ) 

1.392  -5.541-4) 

1.395  -5.581-4) 

1.398  -5 . 62 ( -4 ) 

2.0000 

1.384  -0.00126 

1.386  -0.00127 

1.389  -0.00128 

1.391  -0.00129 

2.2500 

1.370  -0.00180 

1.372  -0.00181 

1.375  -0.00183 

1.377  -0.00184 

2.5000 

1.344  -0.00376 

1.346  -0.00379 

1.348  -0.00381 

1.351  -0.00384 

2.7000 

1.303  -0.00570 

1.305  -0.00574 

1.307  -0.00578 

1.309  -0.00582 

3.0000 

1.293  -0.0955 

1.292  -0.0945 

1.292  -0.0936 

1.291  -0.0926 

3.2000 

1.311  -0.135 

1.309  -0.136 

1.307  -0.137 

1.304  -0.138 

3.3923 

1.352  -0.159 

1.354  -0.161 

1.355  -0.163 

1.357  -0.166 

3.5000 

1.376  -0.158 

1.377  -0.160 

1.378  -0.163 

1.378  -0.165 

3.7500 

1.396  -0.131 

1.399  -0.135 

1.402  -0.138 

1.404  -0.142 

4.0000 

1.398  -0.126 

1.402  -0.129 

1.407  -0.131 

1.411  -0.134 

4.5000 

1.385  -0.120 

1.390  -0.124 

1.395  -0.128 

1.400  -0.132 

5.0000 

1.360  -0.121 

1.366  -0.123 

1.372  -0.125 

1.378  -0.126 

5.5000 

1.337  -0.183 

1.342  -0.184 

1.346  -0.186 

1.351  -0.187 

6.0000 

1.425  -0.195 

1.432  -0.195 

1.439  -0.195 

1.446  -0.195 

6.2000 

1.424  -0.165 

1.431  -0.163 

1.438  -0.161 

1.446  -0.158 

6.5000 

1.370  -0.128 

1.376  -0.122 

1.382  -0.116 

1.388  -0.110 

7.2000 

1.210  -0.176 

1.213  -0.175 

1.216  -0.174 

1.219  -0.172 

7.9000 

1.140  -0.488 

1.141  -0.485 

1.142  -0.482 

1.143  -0.478 

8.2000 

1.200  -0.645 

1.202  -0.652 

1.205  -0.660 

1.207  -0.667 

8.5000 

1.370  -0.755 

1.380  -0.765 

1.390  -0.774 

1.401  -0.784 

8.7000 

1.530  -0.772 

1.556  -0.777 

1.582  -0.782 

1.609  -0.787 

9.0000 

1.650  -0.633 

1.664  -0.630 

1.678  -0.626 

1.691  -0.623 

9.2000 

1.600  -0.586 

1.612  -0.583 

1.624  -0.579 

1.637  -0.576 

9.5000 

1.670  -0.750 

1.684  -0.748 

1.698  -0.747 

1.713  -0.745 

9.8000 

1.910  -0.680 

1.936  -0.650 

1.961  -0.622 

1.987  -0.595 

10.0000 

1.890  -0.455 

1.914  -0.434 

1.938  -0.414 

1.962  -0.395 

10.5910 

1.720  -0.340 

1.730  -0.321 

1.740  -0.304 

1.750  -0.287 

11.0000 

1.670  -0.485 

1.688  -0.479 

1.706  -0.472 

1.724  -0.466 

11.5000 

1.890  -0.374 

1.906  -0.347 

1.922  -0.322 

1.937  -0.299 

12.5000 

1.740  -0.198 

1.748  -0.185 

1.756  -0.172 

1.764  -0.161 

13.0000 

1.690  -0.195 

1.702  -0.183 

1.714  -0.171 

1.725  -0.161 

14.0000 

1.640  -0.195 

1.652  -0.187 

1.663  -0.178 

1.675  -0.171 

14.8000 

1.610  -0.205 

1.616  -0.200 

1.622  -0.194 

1.629  -0.189 

15.0000 

1.590  -0.211 

1.598  -0.207 

1.606  -0.204 

1.613  -0.200 

16.4000 

1.520  -0.414 

1.540  -0.419 

1.560  -0.425 

1.581  -0.431 

17.2000 

1.724  -0.590 

1.757  -0.593 

1.790  -0.596 

1.822  -0.599 

18.0000 

1.950  -0.410 

1.980  -0.369 

2.010  -0.331 

2.041  -0.298 

18.5000 

1.927  -0.302 

1.940  -0.263 

1.954  -0.229 

1.967  -0.200 
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LAMBDA 

(U*> 

300  K 

“r  ni 

290  K 

nr  ni 

280  K 

nr  ni 

270  X 

nr  ni 

20.0000 

1.810  -0.230 

1.820  -0.195 

1.830  -0.165 

1.840  -0.140 

21.3000 

1.790  -0.250 

1.792  -0.235 

1.794  -0.220 

1.797  -0.207 

22.5000 

1.820  -0.290 

1.834  -0.275 

1.849  -0.262 

1.863  -0.248 

25.0000 

1.840  -0.240 

1.858  -0.186 

1.876  -0.144 

1.894  -0.112 

27.9000 

1.780  -0.250 

1.796  -0.188 

1.812  -0.141 

1.828  -0.106 

30.0000 

1.730  -0.290 

1.746  -0.218 

1.762  -0.164 

1.778  -0.124 

35.0000 

1.720  -0.520 

1.736  -0.374 

1.752  -0.269 

1.768  -0.193 

40.0000 

1.940  -0.630 

1.958  -0.524 

1.976  -0.435 

1.994  -0.362 

50.0000 

2.010  -0.650 

2.028  -0.578 

2.046  -0.513 

2.064  -0.456 

60.0000 

2.010  -0.730 

2.028  -0.647 

2.046  -0.574 

2.064  -0.509 

eo.oooo 

2.020  -0.720 

2.036  -0.643 

2.052  -0.575 

2.068  -0.514 

100.0000 

2.020  -0.650 

2.038  -0.584 

2.056  -0.524 

2.074  -0.471 

150.0000 

2.030  -0.560 

2.048  -0.488 

2.066  -0.424 

2.084  -0.369 

200.0000 

2.040  -0.490 

2.056  -0.410 

2.072  -0.342 

2.088  -0.286 

300.0000 

2.060  -0.390 

2.074  -0.297 

2.088  -0.226 

2.102  -0.172 

Table  A-l.  (Cont.) 


i  — “ — - : - : — - - - 1 

LAMBDA 

260  X 

250  X 

240  X 

230  X 

(MM) 

nr  "i 

nr  ni 

nr  ni 

nr  ni 

0.2000 

1.511  -1 . 03 ( -8 ) 

1.515  -1 . 04 ( -8 ) 

1.518  -1.051-8) 

1.521  -1.061-8) 

0.2500 

1.497  -1 . 0 3 ( -8 ) 

1.500  -1 . 04 ( -8  ) 

1.503  -1 . 05 ( -8 ) 

1.507  -1 . 06 ( -8 ) 

0.3000 

1.481  -1 . 03 ( -8 ) 

1.484  -1 . 04  ( -8  ) 

1.488  -1 . 05 ( -8  ) 

1.491  -1.06(-8) 

0.3371 

1.471  -1 . 03 ( -8  ) 

1.474  -1 . 04 ( -8  ) 

1.477  -1 . 05 (-8  ) 

1.480  -1.064-8) 

0.4000 

1.451  -1 . 03 ( -8 ) 

1.454  -1 . 04 ( -8 ) 

1.457  -1 . 05 ( -8 ) 

1.460  -1 . 05 ( -8 ) 

0.4880 

1.443  -1 . 03 ( -8 ) 

1.446  -1.041-8) 

1.449  -1 . 05 ( -8 ) 

1.452  -1 . 05 ( -8  ) 

0.5145 

1.442  -1.031-8) 

1.445  -1 . 04  ( -8  ) 

1.448  -1 . 05 ( -8 ) 

1.451  -1 . 05 ( -8 ) 

0.5500 

1.441  -1.031-8) 

1.444  -1.041-8) 

1.447  -1 . 05 ( -8 ) 

1.450  -1 . 05 ( -8 ) 

0.6328 

1.440  -1 . 52 ( -8 ) 

1.443  -1.531-8) 

1.446  -1 . 54 (-8 ) 

1.448  -1.554-8) 

0.6943 

1.439  -2 . 05 ( - 8  ) 

1.442  -2 . 07 ( -8  ) 

1.445  -2 . 08 ( -8  ) 

1.447  -2.104-8) 

0.8600 

1.436  -1 . 84 ( -7 ) 

1.439  -1.861-7) 

1.442  -1 . 87 ( -7  ) 

1.444  -1.894-7) 

1.0600 

1.431  -1 . 55 ( -6 ) 

1.433  -1.561-6) 

1.436  -1 . 57 ( -6 ) 

1.439  -1.584-6) 

1.3000 

1.421  -1 . 03 ( -5 ) 

1.423  -1.041-5) 

1.426  -1 . 0  5 ( -5 ) 

1.428  -1 .054-5) 

1.5360 

1.413  -1.4H-4) 

1.416  -1.421-4) 

1.419  -1 . 4 3 ( -4 ) 

1.421  -1.444-4) 

1.8000 

1.400  -5 . 66 ( -4 ) 

1.402  -5.701-4) 

1.405  -5 .75 ( -4  ) 

1.407  -5.794-4) 

2.0000 

1.394  -0.00130 

1.396  -0.00131 

1.399  -0.00132 

1.401  -0.00133 

2 .2500 

1.379  -0.00185 

1.382  -0.00187 

1.384  -0.00188 

1.386  -0.00189 

2.5000 

1.353  -0.00387 

1.355  -0.00389 

1.357  -0.00392 

1.359  -0.00395 

2.7000 

1.310  -0.00586 

1.312  -0.00589 

1.314  -0.00594 

1.316  -0.00597 

3.0000 

1.291  -0.0916 

1.290  -0.0907 

1.289  -0.0898 

1.289  -0.0888 

3.2000 

1.302  -0.139 

1.300  -0.140 

1.298  -0.141 

1.296  -0.142 

3.3923 

1.358  -0.168 

1.360  -0.170 

1.362  -0.172 

1.363  -0.175 

3.5000 

1.379  -0.168 

1.380  -0.170 

1.381  -0.173 

1.382  -0.175 

3.7500 

1.407  -0.146 

1.410  -0.150 

1.413  -0.154 

1.416  -0.158 

4.0000 

1.416  -0.137 

1.420  -0.140 

1.424  -0.143 

1.429  -0.146 

4.5000 

1.405  -0.136 

1.410  -0.140 

1.415  -0.144 

1.420  -0.149 

5.0000 

1.384  -0.128 

1.390  -0.130 

1.396  -0.132 

1.402  -0.134 

5.5000 

1.355  -0.189 

1.360  -0.190 

1.365  -0.191 

1.369  -0.193 

6.0000 

1.453  -0.195 

1.460  -0.195 

1.467  -0.195 

1.474  -0.195 

6.2000 

1.453  -0.156 

1.460  -0.154 

1.467  -0.152 

1.474  -0.150 

6.5000 

1.394  -0.105 

1.400  -0.100 

1.406  -0.0952 

1.412  -0.0906 

7.2000 

1.222  -0.171 

1.225  -0.170 

1.228  -0.169 

1.231  -0.168 

7.9000 

1.144  -0.475 

1.145  -0.472 

1.146  -0.469 

1.147  -0.466 

8.2000 

1.210  -0.674 

1.212  -0.682 

1.214  -0.690 

1.217  -0.697 

8.5000 

1.411  -0.794 

1.421  -0.804 

1.431  -0.814 

1.441  -0.824 

8.7000 

1.635  -0.793 

1.661  -0.798 

1.687  -0.803 

1.713  -0.809 

9.0000 

1.705  -0.619 

1.719  -0.616 

1.733  -0.613 

1.747  -0.609 

9.2000 

1.649  -0.572 

1.661  -0.569 

1.673  -0.566 

1.685  -0.562 

9.5000 

1.727  -0.744 

1.741  -0.742 

1.755  -0.740 

1.769  -0.739 

9.8000 

2.012  -0.569 

2.038  -0.544 

2.064  -0.520 

2.089  -0.498 

10.0000 

1.986  -0.377 

2.010  -0.360 

2.034  -0.344 

2.058  -0.328 

10.5910 

1.760  -0.271 

1.770  -0.256 

1.780  -0.242 

1.790  -0.229 

11.0000 

1.742  -0.460 

1.760  -0.454 

1.778  -0.448 

1.796  -0.442 

11.5000 

1.953  -0.278 

1.969  -0.258 

1.985  -0.240 

2.001  -0.222 

12.5000 

1.772  -0.150 

1.780  -0.140 

1.788  -0.131 

1.796  -0.122 

13.0000 

1.737  -0.150 

1.749  -0.141 

1.761  -0.132 

1.773  -0.124 

14.0000 

1.686  -0.163 

1.698  -0.156 

1.710  -0.149 

1.721  -0.143 

14.8000 

1.635  -0.184 

1.641  -0.179 

1.647  -0.174 

1.653  -0.170 

15.0000 

1.621  -0.196 

1.629  -0.193 

1.637  -0.190 

1.645  -0.186 

16.4000 

1.601  -0.436 

1.621  -0.442 

1.641  -0.448 

1.661  -0.454 

17.2000 

1.855  -0.602 

1.888  -0.605 

1.921  -0.608 

1.954  -0.611 

18.0000 

2.071  -0.268 

2.101  -0.241 

2.131  -0.217 

2.161  -0.195 

18.5000 

1.981  -0.174 

1.994  -0.152 

2.007  -0.132 

2.021  -0.115 

Table  A-l .  ( Cont . ) 


: — - - - - - 1 

LAMBDA 

O m) 

260  K 
nr  n. 

2S0  K 

nr  ni 

240  X 
nr  ni 

230  X 

nr  nl 

20.0000 

1.850  -0.118 

1.860  -0.100 

1.870  -0.0847 

1.880  -0.0717 

21.3000 

1.799  -0.194 

1.801  -0.182 

1.803  -0.171 

1.805  -0.160 

22 . 5000 

1.878  -0.236 

1.892  -0.224 

1.906  -0.213 

1.921  -0.202 

25.0000 

1.912  -0.0865 

1.930  -0.0670 

1.948  -0.0519 

1.966  -0.0402 

27.9000 

1.844  -0.0798 

1.860  -0.0600 

1.876  -0.0451 

1.892  -0.0339 

30.0000 

1.794  -0.0930 

1.810  -0.0700 

1.826  -0.0527 

1.842  -0.0396 

35.0000 

1.784  -0.139 

1.800  -0.100 

1.816  -0.0719 

1.832  -0.0517 

40.0000 

2.012  -0.30 

2.030  -0.250 

2.048  -0.208 

2.066  -0.173 

50.0000 

2.082  -0.405 

2.100  -0.360 

2.118  -0.320 

2.136  -0.284 

60.0000 

2.082  -0.451 

2.100  -0.400 

2.118  -0.355 

2.136  -0.314 

80.0000 

2.084  -0.459 

2.100  -0.410 

2.116  -0.366 

2.132  -0.327 

100.0000 

2.092  -0.423 

2.110  -0.380 

2.128  -0.341 

2.146  -0.307 

150.0000 

2.102  -0.322 

2.120  -0.280 

2.138  -0.244 

2.156  -0.212 

200.0000 

2.104  -0.239 

2.120  -0.200 

2.136  -0.167 

2.152  -0.140 

300.0000 

2.116  -0.131 

2.130  -0.100 

2.144  -0.0762 

2.158  -0.0580 

LAMBDA 

220  X 

210  K 

200  X 

190  X 

(«■) 

nr  n. 

nr  ni 

nr  ni 

nr  nA 

0.2000 

0.2500 

0.3000 

0.3371 

0.4000 

0.4880 

0.5145 

0.5500 

0.6328 

0.6943 

0.8600 

1.0600 

1.3000 

1.5360 

1.8000 

2.0000 

2.2500 

2.5000 

2.7000 
3.0000 

3.2000 
3.3923 

3.5000 
3.7500 
4.0000 

4.5000 
5.0000 

5.5000 

6.0000 

6.2000 

6.5000 

7.2000 
7.9000 

8.2000 

8.5000 

8.7000 
9.0000 
9.2000 

9.5000 
9.8000 

10.0000 

10.5910 

11.0000 

11.5000 

12.5000 
13.0000 
14.0000 
14.8000 
15.0000 
16.4000 
17.2000 
18.0000 

18.5000 


1.525 

1.510 

1.494 

1.483 

1.463 

1.455 

1.453 

1.452 

1.451 
1.450 
1.447 
1.442 
1.431 

1.424 
1.410 
1.404 
1.389 
1.361 
1.318 
1.288 
1.293 
1.365 
1.382 
1.418 
1.433 

1.425 
1.408 
1.374 

1.481 

1.482 
1.418 
1.234 
1.148 
1.219 

1.452 
1.740 
1.760 
1.698 
1.784 
2.115 
2.082 
1.800 
1.814 
2.016 
1.804 
1.784 
1.733 
1.660 
1.652 
1.682 
1.986 
2.192 
2.034 


-1 . 07 ( -8 ) 
-1 . 07  { -8 ) 
-1 . 06 ( -8  ) 
-1 . 06 ( -8  ) 
-1 . 06 ( -8 ) 
-1.061-8) 
-1 . 06 ( -8 ) 
-1 . 06  < -8  ) 
-1 . 56 ( -8  ) 
-2.1K-BI 
-1 . 90 ( -7 ) 
-1 . 59 ( -6  ) 
-1 . 06 ( -5 ) 
-1 . 45 ( -4  ) 
-5 . 83 ( -4 ) 
-0.00134 
-0.00191 
-0.00398 
-0.00602 
-0.0879 
-0.143 
-0.177 
-0.178 
-0.163 
-0.149 
-0.154 
-0.136 
-0.194 
-0.195 
-0.148 
-0.0862 
-0.166 
-0.463 
-0.705 
-0.835 
-0.814 
-0.606 
-0.559 
-0.737 
-0.476 
-0.313 
-0.216 
-0.436 
-0.206 
-0.114 
-0.116 
-0.137 
-0.165 
-0.183 
-0.460 
-0.614 
-0.175 
-0.101 


1.528 

1.513 

1.497 

1.486 

1.466 

1.457 

1.456 

1.455 

1.454 

1.453 

1.450 

1.444 

1.434 

1.426 

1.412 

1.406 

1.391 

1.363 

1.320 

1.288 

1.291 

1.366 

1.383 

1.421 

1.438 

1.430 

1.414 

1.378 

1.488 

1.489 
1.424 
1.237 
1.149 
1.222 
1.462 
1.766 
1.774 
1.710 
1.798 
2.140 
2.106 
1.810 
1.832 
2.032 
1.812 
1.796 
1.744 
1.666 
1.660 
1.702 
2.019 
2.222 
2.048 


-1 . 07 ( -8 ) 
-1 . 07 ( -8 ) 
-1 . 07 ( -8) 
-1.071-8) 
-1 . 07 ( -8 ) 
-1 . 07 ( -8  ) 
-1.071-8) 
-1 . 07 ( -8  ) 
-1 . 57 ( -8 ) 
-2.131-8) 
-1.9K-7) 
-1.601-6) 
-1 . 07 ( -5 ) 
-1 . 46  ( -4 ) 
-5 . 87 ( -4 ) 
-0.00134 
-0.00192 
-0.00400 
-0.00605 
-0.0870 
-0.144 
-0.179 
-0.180 
-0.167 
-0.152 
-0.158 
-0.138 
-0.196 
-0.195 
-0.146 
-0.0821 
-0.165 
-0.460 
-0.713 
-0.845 
-0.819 
-0.603 
-0.556 
-0.736 
-0.455 
-0.298 
-0.204 
-0.431 
-0.192 
-0.106 
-0.109 
-0.131 
-0.161 
-0.180 
-0.466 
-0.617 
-0.158 
-0.0878 


1.531 

1.516 

1.500 

1.489 

1.469 

1.460 

1.459 

1.458 

1.457 

1.456 

1.453 

1.447 

1.436 

1.429 
1.415 
1.409 
1.393 
1.366 
1.322 
1.287 
1.289 
1.368 
1.384 
1.424 
1.442 
1.435 
1.420 
1.383 

1.495 

1.496 

1.430 
1.240 
1.150 
1.224 
1.472 
1.792 
1.788 
1.722 
1.812 
2.166 
2.130 
1.820 
1.850 
2.048 
1.820 
1.808 
1.756 
1.672 
1.668 
1.722 
2.052 
2.252 
2.061 


-i . 08 ( -8 ) 
-1 . 08 (-8 ) 
-1.081-8) 
-1 . 08 ( -8  ) 
-1.081-8) 
-1 .08 (-8 ) 
-1 , 08 ( -8 ) 
-1.081-8) 
-1 . 58 ( -8  ) 
-2 . 1 4 ( -8  ) 
-1.931-7) 
-1 . 62 ( -6 ) 
-1.081-5) 
-1.471-4) 
-5.911-4) 
-0.00135 
-0.00193 
-0.00403 
-0.00609 
-0.0861 
-0.145 
-0.182 
-0.183 
-0.172 
-0.156 
-0.163 
-0.140 
-0.197 
-0.195 
-0.144 
-0.0781 
-0.164 
-0.457 
-0.721 
-0.856 
-0.825 
-0.599 
-0.552 
-0.734 
-0.435 
-0.285 
-0.193 
-0.425 
-0.178 
-0.0990 
-0.102 
-0.125 
-0.156 
-0.177 
-0.472 
-0.620 
-0.142 
-0.0765 


1.535 

1.520 

1.503 

1.492 

1.472 

1.463 

1.462 

1.461 

1.460 

1.459 

1.455 

1.450 

1.439 
1.432 
1.417 
1.411 
1.396 
1.368 
1.324 
1.286 
1.287 
1.370 
1.385 
1.427 
1.446 

1.440 
1.426 
1.388 

1.502 

1.503 
1.436 
1.243 
1.151 
1.226 
1.482 
1.818 
1.802 
1.734 
1.826 
2.192 
2.154 
1.830 
1.868 
2.064 
1.828 
1.820 
1.768 
1.678 
1.676 
1.742 
2.085 
2.282 
2.074 


-1 . 09 ( -8 ) 
-1 . 09 ( -8 ) 
-1 . 09 ( -8 ) 
-1 . 09 ( -8  ) 
-1 . 09 ( -8 ) 
-1.091-8) 
-1.091-8) 
-1.091-8) 
-1.601-8) 
-2.161-8) 
-1.941-7) 
-1.631-6) 
-1.081-5) 
-1.481-4) 
-5.951-4) 
-0.00136 
-0.00195 
-0.00406 
-0.00613 
-0.0853 
-0.146 
-0.184 
-0.186 
-0.176 
-0.159 
-0.168 
-0.142 
-0.199 
-0.195 
-0.142 
-0.0744 
-0.163 
-0.453 
-0.729 
-0.867 
-0.830 
-0.596 
-0.549 
-0.733 
-0.416 
-0.272 
-0.182 
-0.419 
-0.165 
-0.0924 
-0.0956 
-0.120 
-0.152 
-0.173 
-0.478 
-0.624 
-0.127 
-0.0667 
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Table  A-l .  ( Cont . ) 


i - 

LAMBDA 

(  L«) 

220  X 
nr  i»i 

210  K 

nr  ni 

200  X 

nr  ni 

190  X 

nr  ni 

20.0000 

1.890  -0.0607 

1.900  -0.0514 

1.910  -0.0435 

1.920  -0.0368 

21.3000 

1.808  -0.150 

1.810  -0.141 

1.812  -0.132 

1.814  -0.124 

22.5000 

1.935  -0.192 

1.950  -0.182 

1.964  -0.173 

1.978  -0.164 

25.0000 

1.984  -0.0312 

2.002  -0.0241 

2.020  -0.0187 

2.038  -0.0145 

27.9000 

1.908  -0.0255 

1.924  -0.0192 

1.940  -0.0144 

1.956  -0.0108 

30.0000 

1.858  -0.0298 

1.874  -0.0225 

1.890  -0.0169 

1.906  -0.0127 

35.0000 

1.848  -0.0372 

1.864  -0.0267 

1.880  -0.0192 

1.896  -0.0138 

40.0000 

2.084  -0.144 

2.102  -0.119 

2.120  -0.0992 

2.138  -0.0825 

50.0000 

2.154  -0.253 

2.172  -0.224 

2.190  -0.199 

2.208  -0.177 

60.0000 

2.154  -0.279 

2.172  -0.247 

2.190  -0.219 

2.208  -0.194 

80.0000 

2.148  -0.292 

2.164  -0.261 

2.180  -0.233 

2.196  -0.209 

100.0000 

2.164  -0.275 

2.182  -0.247 

2.200  -0.222 

2.218  -0.200 

150.0000 

2.174  -0.185 

2.192  -0.161 

2.210  -0.140 

2.228  -0.122 

200.0000 

2.168  -0.117 

2.184  -0.0977 

2.200  -0.0816 

2.216  -0.0682 

300.0000 

2.172  -0.0442 

2.186  -0.0337 

2.200  -0.0256 

2.214  -0.0195 

Table  A-l .  (Cont. ) 


LAMBDA 

( l») 

180  X 

“r  ni 

170  X 

nr  ni 

160  X 

nr  ni 

150  X 

nr  ni 

0.2000 

1.538  -1 . 1 0 ( -8  ) 

1.542  -1.1K-8) 

1.545  -1 . 12 ( -8 ) 

1.549  -1.131-8) 

0.2500 

1.523  -1 . 10 ( -8  ) 

1.526  -l.ll(-B) 

1.529  -1 . 1 2 ( -8 ) 

1.533  -1 . 1 2 ( — 8 ) 

0.3000 

1.506  -1 . 1 0 ( -8 ) 

1.510  -l.ll<-8! 

1.513  -1.111-8) 

1.516  -1 . 12 ( -8 ) 

0.3371 

1.495  -1 . 10 ( -8  ) 

1.499  -1 . 10 ( -8  ) 

1.502  -1.111-8) 

1.505  -1 . 12 { -8 ) 

0.4000 

1.475  -1 . 09 ( -8 ) 

1.478  -1 . 10 { -8 ) 

1.480  -1.111-8) 

1.483  -1 . 12 ( -8 ) 

0.4880 

1.466  -1 . 09 ( -8  ) 

1.469  -1 . 10 ( -8 ) 

1.472  -1.1K-8) 

1.474  -1.121-8) 

0.5145 

1.465  -1 . 09 ( -8 ) 

1.468  -1 . 1 0 ( -8 ) 

1.470  -1.111-8) 

1.473  -1.121-8) 

0 . 5500 

1.464  -1 . 09 ( -8 ) 

1.467  -1.10(-e) 

1.469  -1.111-8) 

1.472  -1.121-8) 

0.6328 

1.463  -1.6K-8) 

1.465  -1 . 62 ( -8 ) 

1.468  -1 . 63 (-8 ) 

1.471  -1.641-8) 

0.6943 

1.461  -2 . 18 ( -8 ) 

1.464  -2.191-8) 

1.467  -2.2K-8) 

1.470  -2.221-8) 

0.8600 

1.458  -1 . 96 ( -7 1 

1.461  -1 . 97  < -7 ) 

1.464  -1 . 99 ( -7 ) 

1.467  -2.001-7) 

1.0600 

1.453  -1.641-6) 

1.455  -1 . 65 ( -6 ) 

1.458  -1 . 66 ( -6 ) 

1.461  -1.681-6) 

1.3000 

1.442  -1.091-5) 

1.444  -1 . 10 ( -5 ) 

1.447  -1.1H-5) 

1.450  -1.121-5) 

1.5360 

1.434  -1.501-4) 

1.437  -1.5K-4) 

1.439  -1 . 52 ( -4 ) 

1.442  -1.531-4) 

1.8000 

1.420  -6 . 00 ( -4 ) 

1.423  -6 . 04  < -4 ) 

1.425  -6 . 08 ( -4  ) 

1.428  -6.121-4) 

2.0000 

1.413  -0.00137 

1.416  -0.00138 

1.418  -0.00139 

1.421  -0.00140 

2.2500 

1.398  -0.00196 

1.401  -0.00197 

1.403  -0.00199 

1.405  -0.00200 

2.5000 

1.370  -0.00408 

1.372  -0.00411 

1.374  -0.00414 

1.376  -0.00417 

2.7000 

1.325  -0.00617 

1.327  -0.00621 

1.329  -0.00625 

1.331  -0.00630 

3.0000 

1.286  -0.0844 

1.285  -0.0835 

1.285  -0.0827 

1.284  -0.0818 

3.2000 

1.285  -0.147 

1.282  -0.148 

1.280  -0.149 

1.278  -0.151 

3.3923 

1.371  -0.187 

1.373  -0.189 

1.374  -0.192 

1.376  -0.194 

3.5000 

1.386  -0.188 

1.386  -0.191 

1.387  -0.194 

1.388  -0.197 

3.7500 

1.430  -0.181 

1.432  -0.186 

1.435  -0.191 

1.438  -0.197 

4.0000 

1.451  -0.162 

1.455  -0.166 

1.460  -0.169 

1.464  -0.173 

4.5000 

1.445  -0.174 

1.450  -0.179 

1.455  -0.185 

1.460  -0.191 

5.0000 

1.432  -0.144 

1.438  -0.146 

1.444  -0.148 

1.450  -0.150 

5.5000 

1.392  -0.200 

1.397  -0.202 

1.401  -0.203 

1.406  -0.205 

6.0000 

1.509  -0.195 

1.516  -0.195 

1.523  -0.195 

1.530  -0.195 

6.2000 

1.510  -0.140 

1.518  -0.138 

1.525  -0.136 

1.532  -0.134 

6.5000 

1.442  -0.0708 

1.448  -0.0674 

1.454  -0.0641 

1.460  -0.0610 

7.2000 

1.246  -0.162 

1.249  -0.161 

1.252  -0.160 

1.255  -0.159 

7.9000 

1.152  -0.450 

1.153  -0.447 

1.154  -0.445 

1.155  -0.442 

8.2000 

1.229  -0.737 

1.231  -0.746 

1.234  -0.754 

1.236  -0.762 

8.5000 

1.492  -0.878 

1.503  -0.889 

1.513  -0.900 

1.523  -0.912 

8.7000 

1.844  -0.836 

1.871  -0.841 

1.897  -0.847 

1.923  -0.853 

9.0000 

1.816  -0.593 

1.829  -0.590 

1.843  -0.587 

1.857  -0.583 

9.2000 

1.746  -0.546 

1.759  -0.543 

1.771  -0.540 

1.783  -0.536 

9.5000 

1.840  -0.731 

1.855  -0.729 

1.869  -0.728 

1.883  -0.726 

9.8000 

2.217  -0.398 

2.243  -0.381 

2.268  -0.364 

2.294  -0.348 

10.0000 

2.178  -0.259 

2.202  -0.247 

2.226  -0.236 

2.250  -0.225 

10.5910 

1.840  -0.172 

1.850  -0.163 

1.860  -0.154 

1.870  -0.145 

11.0000 

1.886  -0.414 

1.904  -0.408 

1.922  -0.403 

1.940  -0.398 

11.5000 

2.080  -0.153 

2.095  -0.142 

2.111  -0.132 

2.127  -0.123 

12.5000 

1.836  -0.0862 

1.844  -0.0804 

1.852  -0.0750 

1.860  -0.0700 

13.0000 

1.832  -0.0896 

1.843  -0.0839 

1.855  -0.0787 

1.867  -0.0737 

14.0000 

1.779  -0.114 

1.791  -0.109 

1.802  -0.105 

1.814  -0.100 

14.8000 

1.684  -0.148 

1.691  -0.144 

1.697  -0.140 

1.703  -0.136 

15.0000 

1.684  -0.170 

1.691  -0.167 

1.699  -0.164 

1.707  -0.161 

16.4000 

1.762  -0.484 

1.783  -0.491 

1.803  -0.497 

1.823  -0.504 

17.2000 

2.118  -0.627 

2.150  -0.630 

2.183  -0.633 

2.216  -0.636 

18.0000 

2.312  -0.115 

2.343  -0.103 

2.373  -0.0926 

2.403  -0.0833 

18.5000 

2.088  -0.0581 

2.101  -0.0507 

2.115  -0.0442 

2.128  -0.0385 
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Table  A-l .  ( Cont . ) 


LAMBDA 

(PM) 

180  X 
nr  ni 

170  X 

nr  ni 

160  X 
nr  n£ 

150  X 

nr  ni 

20.0000 

1.930  -0.312 

1.940  -0.0264 

1.950  -0.0223 

1.960  -0.0189 

21.3000 

1.816  -0.117 

1.819  -0.110 

1.821  -0.103 

1.823  -0.0965 

22.5000 

1.993  -0.156 

2.007  -0.148 

2.022  -0.141 

2.036  -0.134 

25.0000 

2.056  -0.0112 

2.074  -0.00870 

2.092  -0.00674 

2.110  -0.00522 

27.9000 

1.972  -0.00814 

1.988  -0.00612 

2.004  -0.00460 

2.020  -0.00346 

30.0000 

1.922  -0.00957 

1.938  -0.00720 

1.954  -0.00542 

1.970  -0.00408 

35.0000 

1.912  -0.00994 

1.928  -0.00715 

1.944  -0.00514 

1.960  -0.00370 

40.0000 

2.156  -0.0685 

2.174  -0.0570 

2.192  -0.0474 

2.210  -0.0394 

50.0000 

2.226  -0.157 

2.244  -0.140 

2.262  -0.124 

2.280  -0.110 

60.0000 

2.226  -0.172 

2.244  -0.153 

2.262  -0.135 

2.280  -0.120 

80.0000 

2.212  -0.186 

2.228  -0.167 

2.244  -0.149 

2.260  -0.133 

100.0000 

2.236  -0.179 

2.254  -0.161 

2.272  -0.145 

2.290  -0.130 

150.0000 

2.246  -0.106 

2.264  -0.0924 

2.282  -0.0804 

2.300  -0.0700 

200.0000 

2.232  -0.0570 

2.248  -0.0477 

2.264  -0.0399 

2.280  -0.0333 

300.0000 

2.228  -0.0149 

2.242  -0.0113 

2.256  -0.00863 

2.270  -0.00657 

t 


Table  A-2.  Radiative  Parameters  for  the  Proposed  Background 
Stratospheric  Aerosol  Model  at  180  K 


LAMBDA 

( i ia) 

INDEX  OP 
REFRACTION 

Dr  ni 

EXTINCTION 

(km-1) 

SCATTERING 
(km-1 ) 

SINGLE 

ABSORPTION  SCATTERING 
(km-1)  ALBEDO 

ASYMMETRY 

PARAMETER 

0.2000 

1 .5380 

0.0000 

7 . 9699 ( -4 ) 

7.96991-4) 

9.69821-11) 

1.0000 

0.6675 

0.2500 

1.5230 

0 . 0000 

7.76721-4 ) 

7.76721-4) 

7.52431-11) 

1 . 0000 

0.6781 

0.3000 

1.5060 

0.0000 

7.22361-4) 

7.22361-4) 

5.97021-11) 

1.0000 

0.6880 

0 . 3371 

1.4950 

0.0000 

6.70581-4) 

6.70581-4) 

5.17801-11 ) 

1.0000 

0.6926 

0 . 40u0 

1.4750 

0.0000 

5.72431-4 ) 

5.72431-4) 

4.11271-11) 

1.0000 

0.6976 

0.4880 

1.4660 

0.0000 

4.58861-4 ) 

4.58861-4) 

3.17151-11) 

1.0000 

0.6894 

0.5145 

1.4650 

0.0000 

4.29601-4) 

4.29601-4) 

2.95181-11) 

1.0000 

0.6857 

0.5500 

1 . 4640 

0.0000 

3.93521-4) 

3.93521-4) 

2.71001-11) 

1.0000 

0.6799 

0.6328 

1.4630 

0.0000 

3.21881-4) 

3.21881-4) 

3  .  34371-11  ) 

1.0000 

0.6648 

0.6943 

1.4610 

0 . 0000 

2.76591-4) 

2.76591-4  ) 

3.99961-11) 

1.0000 

0.6537 

0.8600 

1  .  4580 

o.onoo 

1.87491-4) 

1.87491-4) 

2.68291-10) 

1.0000 

0.6220 

1.0600 

1 .4530 

0.0000 

1.20171-4 ) 

1.20171-4) 

1.67551-9) 

1 . 0000 

0.5838 

1 . 3000 

1.4420 

0.0000 

7.19021-5) 

7.18921-5) 

8.34281-9) 

1.0000 

0.5391 

1 . 5360 

1.4340 

-0.0002 

4.56321-5) 

4.55401-5) 

9.10621-8) 

0.9980 

o .4yso 

1.8000 

1.4200 

-0.0006 

2.80281-5) 

2.77341-5) 

2.93531-7) 

0.9900 

0.4500 

2 . 0000 

1.4130 

-0.0014 

2.03941-5) 

1.98111-5) 

5.83291-7) 

0.9710 

0.4168 

2.2500 

1.3980 

-0.0020 

1.36661-5) 

1.29481-5) 

7.18211-7) 

0.9470 

0.3778 

2.5000 

1.3700 

-0.0041 

9.36721-6) 

8.05091-6) 

1.31631-6) 

0.8590 

0.3394 

2.7000 

1.3250 

-0.0062 

6.68081-6) 

4.84931-6) 

1.83151-6) 

0.7260 

0.3076 

3.0000 

1.2860 

-0.0844 

2.46621-5) 

2.67051-6) 

2.19911-5) 

0.1080 

0.2600 

3.2000 

1.2850 

-0.1470 

3.77091-5) 

2.39081-6) 

3.53171-5) 

0.0630 

0.2341 

3.3923 

1.3710 

-0.1870 

4.38651-5) 

3.22911-6) 

4.06371-5) 

0.0740 

0.2231 

3.5000 

1.3860 

-0.1880 

4.23651-5) 

3.11621-6) 

3.92491-5) 

0.0740 

0.2160 

3.7500 

1.4300 

-0.1810 

3.73401-5) 

2.93961-6) 

3.44011-5) 

0.0790 

0.2025 

4.0000 

1.4510 

-0.1620 

3.09711-5) 

2.52191-6) 

2.84481-5) 

0.0810 

0.1888 

4.5000 

1 .4450 

-0.1740 

2.85591-5) 

1.66921-6) 

2.68911-5) 

0.0580 

0.1565 

5.0000 

1.4320 

-0.1440 

2,10861-5) 

1.06121-6) 

2.00241-5) 

0.0500 

0.1307 

5.5000 

1.3920 

-0.2000 

2.63171-5) 

6.99191-7) 

2.56171-5) 

0.0270 

0.1066 

6.0000 

1.5090 

-0.1950 

2.19431-5) 

7.54691-7) 

2.11871-5) 

0.0340 

0.0976 

6.2000 

1 . 5100 

-0.1400 

1.53411-5) 

6.34981-7) 

1.47051-5) 

0.0410 

0.0923 

6.5000 

1.4420 

-0.0708 

7.79211-6) 

3.93971-7) 

7.39831-6) 

0.0510 

0.0812 

7 . 2000 

1.2460 

-0.1620 

1.71481-5) 

1.20831-7) 

1.70271-5) 

0.0070 

0.0598 

7.9000 

1.1520 

-0.4500 

4.60081-5) 

2.31771-7) 

4.57761-5) 

0.0050 

0.0460 

8 . 2000 

1 . 2290 

-0.7370 

7.00521-5) 

5.46431-7) 

6.95061-5) 

0.0080 

0.0411 

8.5000 

1  .4920 

-0.8780 

6.42181-5) 

7.03791-7) 

6.35151-5) 

0.0110 

0.0420 

8.7000 

1.8440 

-0.8360 

4.49741-5) 

7.15561-7) 

4.42571-5) 

0.0160 

0.0503 

9.0000 

1.8160 

-0.5930 

3.27141-5) 

4.61271-7) 

3.22521-5) 

0.0140 

0.0499 

9 . 2000 

1.7460 

-0.5460 

3.12251-5) 

3.69891-7) 

3.08571-5) 

0.0120 

0.0463 

9.5000 

1.8400 

-0.7310 

3.64871-5) 

4.46861-7) 

3.60411-5) 

0.0120 

0.0438 

9.8000 

2.2170 

-0.3980 

1.50951-5) 

4.22081-7) 

1.46731-5) 

0.0280 

0.0562 

10.0030 

2.1780 

-0.2590 

1.00921-5) 

3.54741-7) 

9.73771-6) 

0.0350 

0.0534 

10.59. 0 

1.8400 

-0.1720 

8.15151-6) 

1.72481-7) 

7.97891-6) 

0.0210 

0.0388 

11  .ooco 

1.8860 

-0.4140 

1.77101-5) 

1.86811-7) 

1.75231-5) 

0.0110 

0.0360 

11.5000 

2.0800 

-0.1530 

5.54751-6 ) 

1.73801-7) 

5.37381-6) 

0.0310 

0.0378 

12.5000 

1.8360 

-0.0862 

3.46291-6) 

8.50251-8) 

3.37781-6) 

0.0250 

0.0278 

13.0000 

1.8320 

-0.0896 

3.45301-6) 

7.20771-8) 

3.38101-6) 

0.0210 

0.0256 

14.0000 

1.7790 

-0.1140 

4.20231-6) 

4.87901-8) 

4.15351-6) 

0.0120 

0.0214 

14.8000 

1.6840 

-0.1480 

5.51841-6) 

3.26171-8) 

5.48581-6) 

0.0060 

0.0182 

15.0000 

1.6840 

-0.1700 

6.24321-6) 

3.13511-8) 

6.21181-6) 

0.0050 

0.0177 

16.4000 

1.7620 

-0.4840 

1.49211-5) 

3.44901-8) 

1.48871-5) 

0.0020 

0.0149 

17.2000 

2.1180 

-0.6270 

1.34281-5) 

4.51551-8) 

1.33831-5) 

0.0030 

0.0161 

18.0000 

2.3120 

-0.1150 

2.12961-6) 

3.53601-8) 

2.09421-6) 

0.0170 

0.0172 

18.5000 

2.0880 

-0.0581 

1.26191-6) 

2.49831-8) 

1.23691-6) 

0.0200 

0.0145 

20.0000 

1.9300 

-0.3120 

6.93021-6) 

1.63571-8) 

6.91401-6) 

0.0020 

0.0112 

21 . 3000 

1.8160 

-0.1170 

2,70311-6) 

9.63971-9) 

2.69351-6) 

0.0040 

0.0094 
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Table  A-2.  (Continued) 


LAMBDA 

(Wa) 

INDEX  OP 
REFRACTION 

nr  ni 

EXTINCTION 
(k*  ) 

SCATTERING 

Hub"1) 

ABSORPTION 

Oca"1) 

SINGLE 

SCATTERING 

ALBEDO 

ASYMMETRY 

PARAMETER 

22.5000 

1.9930 

-0.1560 

2 . 9443  ( -6  ) 

1.02561-8) 

2.93401-6) 

0.0030 

0.0092 

25.0000 

2.0560 

-0.0112 

1 . 8728 ( -7 ) 

7.11291-9) 

1.80171-7) 

0.0380 

0.0078 

27.9000 

1.9720 

-0.0081 

1 . 2970 ( -7  ) 

4.09421-9 ) 

1.25611-7) 

0.0320 

0.0060 

30.0000 

1.9220 

-0.0096 

1 . 4588 ( -7 ) 

2.84411-9) 

1.43031-7) 

0.0190 

0.0050 

35.0000 

1.9120 

-0 . 0099 

1 . 2976 ( -7 ) 

1.50901-9) 

1.28251-7) 

0.0120 

0.0037 

40.0000 

2.1560 

-0.0685 

6 . 3 1 8 7 ( -7  ) 

1  .21241-9) 

6.30661-7) 

0.0020 

0.0032 

50 . 0000 

2.2260 

-0 .1570 

1 . 0877 ( -6 ) 

5.38861-10) 

1.08721-6) 

0.0000 

0.0021 

60 . 0000 

2.2260 

-0 .1720 

9.91591-7) 

2.60381-10) 

9.91321-7) 

0.0000 

0.0015 

80.0000 

2.2120 

-0.1860 

8 . 1250 ( -7  ) 

8.14701-11) 

8.12421-7) 

0.0000 

0.0008 

100.0000 

2 .2360 

-0.1790 

6 . 1 3  1 7  (  -  7  ) 

3.40561-11) 

6.13141-7) 

0.0000 

0.0005 

150 . 0000 

2.2460 

-0.1060 

2 . 4  06 1  ( - 7  ) 

6.70021-12) 

2.40591-7) 

0.0000 

0.0002 

200 . 0000 

2.2320 

-0.0570 

9 . 8257 ( -8  ) 

2.08131-12) 

9.82551-8) 

0 . 0000 

0.0001 

300.0000 

2.2280 

-0.0149 

1  .  7 1 8  5  (  -  8  ) 

4.08711-13) 

1.71851-8) 

0.0000 

0.0001 

A-ll 


Table  A-3.  Radiative  Parameters  for  the  Proposed  Background 
Stratospheric  Aerosol  Model  at  250  K 


INDEX  OF  SINGLE 


LAMBDA 

( u“) 

REFRACTION 

nr  ni 

EXTINCTION 
(k«-i ) 

SCATTERING 

Oun'1) 

ABSORPTION 

Oun_i) 

SCATTERING 

ALBEDO 

ASYMMETRY 

PARAMETER 

0.2000 

1.5150  0.0000 

8 . 4896 ( -4  ) 

8.48961-4) 

9.77981-11) 

1.0000 

0.6811 

0.2500 

1.5000  0.0000 

8 . 1890 ( -4  ) 

8.18901-4) 

7.35541-11) 

1.0000 

0.6927 

0.3000 

1.4840  0.0000 

7 . 5494  ( -4  ) 

7.54941-4) 

5.92121-11) 

1.0000 

0.7014 

0.3371 

1.4740  0.0000 

6 . 9694  ( -4  ) 

6.96941-4) 

5.15581-11) 

1.0000 

0.7050 

0.4000 

1.4540  0.0000 

5 . 8792 ( -4  ) 

5.87921-4) 

4.05761-11 ) 

1.0000 

0.7093 

0.4860 

1.4460  0.0000 

4.67221-4) 

4.67221-4) 

3.15621-11  ) 

1.0000 

0.6995 

0.5145 

1.4450  0.0000 

4  .  3658  ( -4 ) 

4  .  36581  -4  ) 

2.95441-11) 

1.0000 

0.6952 

0 .5500 

1.4440  0.0000 

3 . 9860 ( -4  ) 

3.98601-4) 

2.70821-11 ) 

1.0000 

0.6892 

0.6328 

1.4430  0.0000 

3.24061-4) 

3.24061-4) 

3.32681-11) 

1.0000 

0.6732 

0.6943 

1.4420  0.0000 

2 . 7826 ( -4  ) 

2.78261-4) 

3.97981-11) 

1 .0000 

0.6611 

0.8600 

1.4390  0.0000 

1 .87261-4  ) 

1  .  8726  1  -4  ) 

2.67861-10) 

1.0000 

0.6279 

1  .0600 

1.4330  0.0000 

1  .  1 87 9 ( -4  ) 

1  .  1879  1  -4  ) 

1 . 68251-9) 

1 .0000 

0.5883 

1.3000 

1.4230  0.0000 

7 . 0870 ( -5 ) 

7.08621-5) 

8.44081-9) 

1.0000 

0.5417 

1  .  5360 

1.4160  -0.0001 

4.49881-5) 

4.48961-5) 

9.17461-8) 

0.9980 

0.4971 

1 .8000 

1.4020  -0.0006 

2.75081-5) 

2.72121-5) 

2.97661-7) 

0.9890 

0.4498 

2.0000 

1.3960  -0.0013 

2 . 0088 ( -5  ) 

1.94911-5) 

5.96481-7) 

0.9700 

0.4160 

2.2500 

1.3820  -0.0019 

1.34911-5) 

1.27571-5) 

7.34081-7) 

0.9460 

0.3765 

2.5000 

1 .3550  -0.0039 

9.27101-6) 

7.92521-6) 

1.34601-6) 

0.8550 

0.3377 

2.7000 

1.3120  -0.0059 

6.65421-6) 

4.77861-6) 

1.87561-6) 

0.7180 

0.3060 

3.0000 

1.2900  -0.0907 

2.81581-5) 

2.94781-6) 

2.52101-5) 

0.1050 

0.2598 

3.2000 

1 .3000  -0.1400 

3.85761-5) 

2.73901-6) 

3.58381-5) 

0.0710 

0.2365 

3.3923 

1 .3600  -0.1700 

4.30021-5) 

3.21161-6) 

3.97901-5) 

0.0750 

0.2235 

3.5000 

1 . 3800  -0 .1700 

4.13281-5) 

3.17221-6) 

3.81561-5) 

0.0770 

0.2170 

3.7500 

1.4100  -0.1500 

3.37121-5) 

2.81301-6) 

3.09001-5) 

0.0830 

0.2028 

4.0000 

1.4200  -0.1400 

2.91061-5) 

2.34181-6) 

2.67661-5) 

0.0800 

0.1868 

4 . 5000 

1.4100  -0.1400 

2.51381-5) 

1.49941-6) 

2.36381-5) 

0 . 0600 

0.1544 

5.0000 

1.3900  -0.1300 

2.07441-5) 

9.38141-7) 

1.98061-5) 

0.0450 

0.1276 

5.5000 

1.3600  -0.1900 

2.71621-5) 

6.46321-7) 

2.65161-5) 

0.0240 

0.1048 

6.0000 

1.4600  -0.1950 

2.40641-5) 

6.88741-7) 

2.33761-5) 

0.0290 

0.0947 

6.2000 

1.4600  -0.1540 

1.84311-5) 

5.79021-7) 

1.78521-5) 

0.0310 

0.0894 

6.5000 

1.4000  -0.1000 

1.18191-5) 

3.60041-7) 

1.14591-5) 

0.0300 

0.0792 

7.2000 

1.2250  -0.1700 

1.94411-5) 

1 .19031-7) 

1.93221-5) 

0.0060 

0.0592 

7.9000 

1.1450  -0.4720 

5.19061-5) 

2.69601-7) 

5.16361-5) 

0.0050 

0.0456 

8.2000 

1.2120  -0.6820 

7.01281-5) 

4.98981-7) 

6.96301-5) 

0.0070 

0.0415 

8.5000 

1.4210  -0.8040 

6.73801-5) 

6.36301-7! 

6.67441-5) 

0.0090 

0.0415 

8.7000 

1.6610  -0.7980 

5.36241-5) 

6.52521-7) 

5.29721-5) 

0.0120 

0.0457 

9.0000 

1.7190  -0.6160 

3.90941-5) 

4.59701-7) 

3.86341-5) 

0.0120 

0.0468 

9.2000 

1.6610  -0.5690 

3.70921-5) 

3.69761-7) 

3.67221-5) 

0.0100 

0.0438 

9.5000 

1.7410  -0.7420 

4.28381-5) 

4.50841-7) 

4.23861-5) 

0.0110 

0.0412 

9.8000 

2.0380  -0.5440 

2.47201-5) 

4.16801-7  ) 

2.43041-5) 

0.0170 

0.0488 

10.0000 

2.0100  -0.3600 

1.67401-5) 

3.28781-7) 

1.64111-5) 

0.0200 

0.0475 

10.5910 

1.7700  -0.2560 

1.35451-5) 

1.71561-7) 

1.33731-5) 

0.0130 

0.0370 

11.0000 

1.7600  -0.4540 

2.28501-5) 

1.76551-7) 

2.26741-5) 

0.0080 

0.0333 

11 . 5000 

1  .  9690  -0.2580 

1.07091-5) 

1.67351-7) 

1.05421-5) 

0.0160 

0.0352 

12.5000 

1.7800  -0.1400 

6.21021-6) 

8 .3484  1-8  ) 

6.12661-6) 

0.0130 

0.0269 

13.0000 

1.7490  -0.1410 

6.13821-6) 

6.71221-8) 

6.07101-6) 

0.0110 

0.0244 

14.0000 

1.6980  -0.1560 

6.52001-6) 

4.51881-8) 

6.47481-6) 

0.0070 

0.0205 

14.8000 

1.6410  -0  1790 

7.36541-6) 

3.23741-8) 

7.33301-6) 

0.0040 

0.0177 

15.0000 

1  .6290  -0.1930 

7.89921-6) 

3.02021-8) 

7.86901-6) 

0.0040 

0.0172 

16 . 4000 

1.6210  -0.4420 

1.64081-5) 

2.87941-8) 

1.63791-5) 

0.0020 

0.0139 

17.2000 

1.8880  -0.6050 

1.68511-5) 

•3.97001-8) 

1.68111-5) 

0.0020 

0.0142 

18.0000 

2.1010  -0.2410 

5.58001-6) 

3.16321-8) 

5.54821-6) 

0.0060 

0.0153 

18.5000 

1.9940  -0.1520 

3.75481-6) 

2.41541-8) 

3.73081-6) 

0.0060 

0.0137 

20.0000 

1.8600  -0.1000 

2.54761-6) 

1.42461-8) 

2.53341-6) 

0.0060 

0.0109 

21.3000 

1.8010  -0.1820 

4.53841-6) 

1.03411-8) 

4.52801-6) 

0.0020 

0.0093 
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INDEX  OP  SINGLB 

LAMBDA  REFRACTION  EXTINCTION  SCATTERING  ABSORPTION  SCATTERING  ASYMMETRY 
(pa)  n  (lea  )  Oca"1)  Oca"1)  ALBEDO  PARAMETER 


22.5000 

1.8920 

-0.2240 

4 .8942 ( -6 ) 

9 . 8082 ( -9 ) 

4.88441-6) 

0.0020 

0.0087 

25.0000 

1.9300 

-0.0670 

1 . 2850  < -6 ) 

6 . 4360 ( -9  ) 

1.27861-6) 

0.0050 

0.0073 

27.9000 

1.8600 

-0.0600 

1 . 089  3 ( -6 ) 

3 . 7050 ( -9 ) 

1.08561-6) 

0.0030 

0.0056 

30.0000 

1 . 8100 

-0 . 0700 

1 . 2287 ( -6  ) 

2 . 5440 ( -9  ) 

1.22611-6) 

0.0020 

0.0047 

35.0000 

1.8000 

-0.1000 

1 . 51 28 ( -6 ) 

1 . 3574  ( -9  > 

1.51151-6) 

0.0010 

0.0035 

40.0000 

2.0300 

-0.2500 

2 . 71 68 ( -6  ) 

1.18001-9) 

2.71561-6) 

0.0000 

0.0030 

50.0000 

2.1000 

-0.3600 

2 . 9246 ( -6  ) 

5.48421-10) 

2.92401-6) 

0.0000 

0.0019 

60.0000 

2.1000 

-0.4000 

2 . 6962 ( -6  ) 

2.70041-10) 

2.69581-6) 

0.0000 

0.0013 

80.0000 

2.1000 

-0.4100 

2 .0692  (-6  ) 

8.58461-11 ) 

2.06921-6) 

0.0000 

0.0008 

100.0000 

2.1100 

-0.3800 

1 . 5255 ( -6 ) 

3.48821-11) 

1  .  5254  1-6  ) 

0.0000 

0.0005 

150.0000 

2 .1200 

-0.2800 

7 . 4886  { -7 ) 

6.63201-12) 

7.48861-7) 

0.0000 

0.0002 

200.0000 

2.1200 

-0.2000 

4.03001-7) 

2.03901-12 ) 

4.03001-7) 

0.0000 

0.0001 

300.0000 

2.1300 

-0.1000 

1 . 3 3 68 ( -7  ) 

3.97941-13) 

1.33681-7) 

0.0000 

0.0001 

*US  GOVERNMENT  PRINT  INC,  OF  F  ICE 
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